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Abstract 

A comprehensive set of experiments was performed to study stimulated Brillouin 
scattering (SBS) in a CO, laser-plasma interaction. The facility consisted of a TEA CO, 
laser producing either 1 GW, 40 ns FWHM gain switched laser pulses or 1.7 GW, 2 ns 
FWHM injection mode locked laser pulses which were focussed to 100 micron spot size 
onto an oxygen gas jet target. The measured average oxygen plasma conditions during 
the SBS were N, =5x 10% cm, Be = 120 eV, Z= 6 with the gain switched laser 
and) nig Sexe Oven? Ve =150 eV, Z=3 with the injection mode locked laser. 
Strong ion wave damping, vehog - 0(1) | was deduced for both cases from 
backscattered light spectra, reflected power and ruby Thomson scattering probing of the 
SBS driven ion modes. Two ion heating models were developed to interpret the 
convective SBS instability growth due to the large linear ion Landau damping. An ion 
trapping and heating description whereby only a fraction of the plasma ions was heated 
to temperature O( ZT, ) was applicable to the 2 ns pump laser while a model 
incorporating energy balance with free streaming ions described the 40 ns pump laser 
under conditions of moderate SBS reflectivity. An 80 micron SBS interaction region 
measured by temporally and spatially resolved ruby Thomson scattering was explained 
through axial plasma density gradients. 

An important experiment demonstrated SBS suppression with a multiline CO, 
pump laser. This improved laser-plasma coupling is highly desirable for direct drive 
laser-fusion designs. The 14% single line laser reflectivity was reduced to less than 2% 
when the gain switched laser oscillated on two or more equal intensity lines with the line 


separation, Aw , greater than the homogeneous plasma growth rate, Y = 107% 


0 
s?. Anomalous SBS mode coupling took place when Aw < Wie and the SBS 


reflectivity was not severely reduced. 
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1. Introduction 


1.1 Motivation for Research into Laser-Plasma Interactions 

Current research in direct drive laser fusion has revealed the complex behaviour 
of laser—plasma interactions. Ideally, when a suitable fuel pellet is irradiated with high 
enough laser light intensity, it will be isentropically compressed and the ions thermally 
heated to thermonuclear burn resulting in large energy yields [1]. However, in both the 
overdense core of the irradiated target and the coronal plasma, many instabilities and 
laser—plasma effects take place to prevent such an event. Of interest here are those 
parametric instabilities which can occur in the ablated coronal plasma surrounding the 
target and consequently degrade the laser—plasma coupling. One instability, the photon to 
plasmon plus phonon parametric decay, has been postulated to be a useful absorptive 
instability to enhance the laser heating of the target. Unfortunately, this is only one of 
four similar resonant normal mode triplets which occur in the underdense plasma of the 
target corona. Stimulated Raman scattering (SRS) of the incident pump laser produces 
reflected photons and Langmuir waves in plasma density < 1/4 no , where We 
is the critical density for the particular laser frequency. Another parametric instability, 
stimulated Brillouin scattering (SBS) couples two photons to the ion acoustic wave. Both 
SRS and SBS are reflective type instabilities where a large fraction of the incident light 
energy is converted to reflected light. At the high laser intensities required to initiate 
fusion, these processes may be severe, even at densities much less than We , where 
the absorptive decay do not occur. SBS is especially recalcitrant as it can occur in any 
density, Urs n= Ce , and it has a low threshold in most laser plasmas. The 
objective of this thesis research was to study stimulated Brillouin backscatter in a CO, 


laser-plasma interaction and experimentally investigate the plasma and laser effects 


which may limit its growth. 
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1.2 Early work in Stimulated Brillouin Scattering 

The first theoretically predicted parametric decay in a plasma was a Langmuir 
wave pumping another langmuir wave and an ion acoustic wave (Oraevsky and Sagdeev 
[2]). It was nine years later when this particular decay was experimentally observed in a 
laboratory plasma (Franklin et al. [3]). Since then, other parametric instabilities including the 
photon decays, SRS, SBS and two plasmon decay have been studied both theoretically 
and experimentally. All these parametric instabilities are normal mode resonant triplets 
where energy exchange takes place with frequency and wavenumber matching of the 
participating waves satisfying their respective plasma dispersion relations. Written 


explicitly, these are: 


= AY 
"0. atuckote 12 
= + 
SG ee 
ex(w;sk;) = 0 = (Oeil 
where © is the wave frequency, k is the wave vector and €; = 9 is the 


‘ 
dispersion relation for a particular wave. The subscripts denote the participating waves 


with O identifying the externally provided pump wave. It is clearly seen by multiplying 
the frequency and wavevector sum rules by Planck’s constant, Nh , that energy and 
momentum are conserved as the pump wave decays to the two product waves. 

Early theoretical work by Nishikawa [4] generalized the parametric process for 
infinite homogeneous plasma. Here, two oscillators were coupled through a third pump 
field and the threshold and growth rates were derived. Perkins and Flick [5] further 
considered the decay of an electromagnetic wave including wavenumber and frequency 
mismatches due to plasma inhomogeneities. Rosenbluth [6] clarified this situation by 
solving the WKB coupled mode equations through Laplace transform methods arriving at 
expressions for absolute growth in an inhomogeneous plasma. Rosenbluth et al. [7] later 
extended this work by discussing the evolution of temporal growth to convective 
saturation. 

Pesme et al. [8] presented a particularly lucid exposition of the coupled equations 


for constant pump which described both absolute and convective wave growth as 
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affected by damping in cases of a homogeneous plasma, finite length plasma or 
inhomogeneous plasma. Their homogeneous bounded plasma equations describing the 


Brillouin instability are representative of the WKB approach: 


* 
ga_/at + c_a_/ax = yn 


* * * 
on /ot + can [o%% aaa ta 


0} =- 
where C_ is the phase velocity of the backscattered electromagnetic wave with 
action a_ , and C,. is the ion sound speed associated with the ion fluctuation 
n . The coupling coefficient, vere is proportional to the pump electromagnetic 


field and ie is the ion wave damping. Upon solving for a dispersion relation where 


Nger 0 , absolute instability in homogeneous plasma was predicted if the 
, 1/2 

plasma dimension exceeded the gain length m(c_c.) / fey, . Wave damping 
imposed the usual instability threshold condition, — y, > rity) is where = T_ 


is the electromagnetic wave damping. Also, because of wave energy convection out from 
the bounded plasma, convective instability took place if 


Wz 


Abe (y./2)(c_/c.) . A result which is directly applicable to many 


6) 

stimulated Brillouin experiments came from ignoring the time derivatives and assuming 

strong ion wave damping. In this case, the instability growth length Cais 

was derived from the coupled equations. Finally, the case of inhomogeneous plasma 

showed that the characteristic dimension of the unstable region was = (k' ye 

where k' is the spatial derivative of the wave number mismatch, 
Ms (ko j 4 x ko) : 

A significant departure from the above descriptions came with Drake et al. [9] 
using the plasma fluid equations to derive a general dispersion relation for parametric 
instabilities. The earlier WKB models required that the spatial and temporal wave growths 
be small compared to the real wavenumber and frequencies. Drake et al. imposed no 
restriction on the pump field amplitude (within nonrelativistic theory) whereby the decay 
waves were no longer necessarily normal modes of a field free plasma and could 


therefore have growth rates comparable to their wavenumbers and frequencies. Forslund 


et al. [10] extended this work with examples of Brillouin and Raman scattering in both 
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fluid and kinetic limits of the plasma susceptibilities. In a second paper, Forslund et al. [1 1] 
showed computer simulation results and discussed nonlinear effects in the instability 
growth and saturation. Later work by Cohen and Max [12] generalized the studies of 
various parametric instabilities, including stimulated Brillouin scattering through a Green's 
function formulation. 

Experimental studies of laser induced plasma parametric instabilities are recent 
when compared to earlier observations of Brillouin and Raman scattering in other 
materials. While Brillouin scattering from phonons in solids was detected by Chiao et al. in 
1964 [13], it was not until the early 1970's that Brillouin scattering measurements in 
plasma were made. Both solid and gas targets have been used in the parametric instability 
research with high power Nd:glass and carbon dioxide lasers. The initial Brillouin 
experiments at low laser intensities and plasma electron densities showed fractions of a 
percent SBS reflectivity [14], but with increased laser power and plasma densities, 
reflectivities approaching unity have been measured [15]. It is clear from these high 
reflectivities that the seriousness of stimulated Brillouin scattering in laser fusion 


experiments entitles it to intensive theoretical and experimenta! research. 


1.3 Objective of this Work 

This thesis attempts to sort out the dominant laser and plasma effects which 
saturated the stimulated Brillouin backscatter and limited its growth rate in a series of 
experiments. In the experiments, a free oxygen jet was irradiated with either a long pulse, 
gain switched CO, laser pulse or a short pulse, injection mode locked (IML) CO, laser 
pulse. The long pulse laser was used to produce a well characterized; highly ionized 
plasma from which the laser light could then scatter. The short pulse laser generated fast 
rising, high intensity light pulses where comparisons to the gain switched laser results 
revealed the importance of different ion heating mechanisms in saturating the Brillouin 
instability. In particular, ion heating is examined and found to explain both the instability's 
temporal and spatial growth. 

Chapter 2 summarizes theoretical considerations of stimulated Brillouin scattering 
in laser-plasma experiments. Both the generalized dispersion relation and nonlinear 


coupled mode equations are examined. The coupled mode approach is shown to be 
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particularly amenable to incorporating a wide variety of problems in both spatial and 
temporal instability growth. With these equations, ion heating and ion trapping models are 
developed. 

Apparatus and diagnostic techniques are fully explained in Chapter 3. This is an 
important chapter describing the IML laser system, laser focussing and the oxygen gas 
target. Also, the infrared light detection equipment used to monitor the Brillouin 
backscattered light and the ruby laser Thomson scattering method are described. 

The experimental results and deductions made from them are contained in 
Chapter 4. The plasma parameters are measured through interferometry and X-ray 
measurements. Next, preliminary observations on the Brillouin instability are given. The 
remainder of the chapter presents measurements of the backscattered light and of the 
ion sound waves of the Brillouin instability. These results were interpreted to show ion 
trapping and heating and plasma density inhomogeneity as the major contributors 
inhibiting the instability's growth. An important experiment was also performed to 
demonstrate Brillouin reflectivity reduction with a multiwavelength pump laser. Concluding 
remarks in Chapter 4 propose several mechanisms which could explain the sudden 


termination of the Brillouin scattering in the laser—plasma interaction. 


lease ei anion o ule soime 
tra nner sie 0! Ss erat 5 oir 


Me aiceia? .£ rgigend) ci cenialgas: yi? as seuss 
em CSpyxo ac? ers triesuse?7 +28F matey: oa 
NUwTS att worry ST 26a THEMGN PS sy 
Soe oeey oe novtign ak la B iaeenantt = 
BRIA K Bre wierowe isin, Cowon oS ueabir ae Pasty 
eT ine oi eta yilirefon mivoiad Aivt pray span 
at Io ons fmol pajyaiecursed eh rt ainamSsr imam | 
Aon wtte.wt HaTsie ei7 sisw Mbess azont Price, 
SUGVULOW » THEI ofl <o Visteoseieiet yi} ote 
ot partiohedq oale asw mernnisdxa inetaggm Abe 
ptubulsro. Iseshiarba 4 fHfgndibvew! rhage a iw Pctiubete 
nbtbue ady nistoxs eves nai etnieinatiaame iajevaa. wang. si 


SitoeveIni sroasig~iseal Orit (h gn yilisae null 


2. Theory 


2.1 Purpose of this Chapter 

In this chapter, the theory of Stimulated Brillouin (Compton) scattering is studied 
with emphasis on those aspects relevant to the experiments discussed in Chapter 4. 
Geometrical effects on the instability are ignored and all theories will be explained 
through one dimensional models. Further details will be given in the analysis of the 
experiments based on the concepts presented here. Introductory remarks on Brillouin 
scattering have already been given in the first chapter and will not be repeated in this 
chapter. 

This chapter is divided into three sections. First, the Brillouin dispersion relation 
derived from Maxwell's equations, Poisson's equation and Viasov’s equation is examined 
in both fluid and kinetic limits for temporal and spatial growth. These results apply to low 
amplitude waves and describe the instability’s strength and bandwidth before nonlinear 
effects modify its behaviour. These new normal modes are not necessarily the 
eigenmodes of the field free plasma. The second section of this chapter is centered on 
examining first order coupled wave equations derived from fluid and Poisson's equations. 
These equations easily describe the weakly coupled temporal growth and _ spatial 
distribution of the participating waves as well as provide a simplified insight into 
distinguishing between convective and absolute wave growth. Nonlinear behaviour in 
Brillouin scattering is examined in this section. The coupled mode equations are solved 
for the simplest case where direct ion heating and pump depletion control the instability’s 
spatial growth. An added benefit of the coupled equations is the straightforward 
extention of this model to the case of a multiwavelength laser pump which will be shown 
to decrease the backscattered light. A second simple model, that of ion trapping, is also 
incorporated into the coupled mode equations and provides an alternate description of 
the spatial growth. All these simplified analyses are phenomenological descriptions rather 
than self-consistent solutions to the wave amplitudes and plasma particle energy. The 
latter approach is generally not tractable except through particle simulations, whereas a 
phenomenological approach allows a reasonable physical picture of the different plasma 


effects. 


by enero Ai Reeeny oie snermigans sit ot 2 eal Shah 
ha as 


lidatari. 
tygirinlcans ed’ liiw asnosm. a ong baron! oi. piliga sau 


| a ee 
ayy 70s Were ani fil Riel ad ane mere: “ene ak ‘writ 
| ‘ 
neti od sh eria Aerabboun $190 pales, aryeotiog | wi 
iar 
aid? ok Sareege7 ed 400 Lin’ pas yeIgeHs toh ‘anit Ai rahe ona 


: a 
Ties? MOlstTemelD atl CLenil ai! apt 307, eng) jose set eh 


ms 


a - bea ; 
2) noOvEuUne. 2 VO8 lly 70F agita pa Bndcatels ahh is vee . 


SSrimexs 
(ro? of Vigds etues) seanT wor In odie. 6) i> roe sah 


sariincd soled rthiwhiase kare Nips: ethan a ing 


au 


issage2sn ton gis ce lammi07 5 Muar, -aaent 4d 
vt or] ok 


iO OSsIsiNes 21 “eras ait 1 nore Wiuoas “SAT. babel 1° ai 


nC 60s @ 7102810 q OGe & ri mot Baniwo ahsireund evewt J 


< 


siege ons wo jstoqnay pelgues: “ylaady "melt, ie 
ae ; , 


‘,' 


otal ingieth: esi iiqenta s abivesg 2) Hav as [eovgw Ms 
ni susivedsd IgpTiinol Atal eg avaiae shat teat 
savin S16 ndisups) abom balques nt atead yt ne é . 
a yiilldara oct loitos fwitaigat amniua bbe Boies | Ae! Tagyib: 236 IW Sabo Tepleirig e 
Viswiciingiase’ Sri) ogi iat et Sf ae at \eiitaned nsbbe oA ~ 
wore So ii risiriw optics ‘baal tan . a ng di iQ? istoaim emt 
TAGs Briagey Hot to tert ietaual 2 bnooee A fro berertasaoad ¢ 
- ¥o notaivseek: siamese na zouivang Bais dhowenps bom lasiauioa, ‘att on 
eras sitcaesa oeoorancnqtn seeyana Daitiqnie: oka iit ow 


2.2 Dispersion Relations 
The usual two component plasma dispersion relation for low amplitude 
electrostatic waves in the absence of other strong electrostatic or electromagnetic 


waves is well known, ie. € = | + el ahs O where xX; and x 


1 e 
are the ion and electron susceptibilities respectively. For a particle velocity distribution 


ur (u) and for a wave with phase velocity w/k 
w 2 of ./ du 


J k2 u - w/k 


where Wis the wave frequency and k is its wavenumber. The free plasma 
eigenmodes are those which are weakly damped. For the ion acoustic wave, the solution 


to this dispersion relation is the familiar 


= = + js ec 2 - 
w/k. Cc. [(ZT,, 37, )/M. J /( a aK An ) EqhiZ 
where Z_ is the ionization state, Ta) the electron (ion) temperature, 
M. the ion mass and hy the plasma Debye length. The imaginary 


frequency component is given by the Landau formula (neglecting electron contribution) 


Im w = w.(n/8)'/2(21./T,)3/@exp(-27,/2T, - 3/2) Eq 23 


7s 

In the presence of a large amplitude electromagnetic wave, strong parametric — 
coupling in the plasma can take place between it and electrostatic or other 
electromagnetic waves. The dispersion relation is modified from above since the solution 
can now admit temporal and spatial growth of the daughter waves at the expense of the 
pump field energy. These parametric instabilities include Raman and Brillouin scattering of 
the incident pump wave, the latter being of greatest interest here. 

The dispersion relation for parametric instabilities has been derived many times 
before (Forslund et al. [10], Drake et al. [9], Cohen and Max [12]) and its derivation will be 
outlined here. The incident linearly polarized electromagnetic pump wave satisfies the 


usual dispersion relation wi = wp” ee od a since the temporal and spatial 
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growth rates of the instability are still small compared to ore and ko 


where %5 is the pump frequency, ko is its free space wavenumber and 


Bs, 2 1/2 
or (ne /e om) 
sideband modes, (w, sk, ) Se Astae ate sk +t ko) generated, together with 


is the plasma frequency. The wave equation for the 


the density fluctuations, én. -(t) « exp(i(k-x - wt)) when Fourier 


transformed yields 


Eq. 2.4 


where EY is the sideband electric field and an is that of the pump. 
This equation describes the forcing function acting on the homogeneous wave equation. 
Equating Eon to O reduces the description to that of a freely propagating 
electromagnetic wave in the plasma. Also note at this point that the resonance conditions, 

Re(w) = Re(w, + w) Re(k.) = Re(k, oaky) apply, making this a 


linear model without higher harmonic wave coupling. From Vlasov and Poisson's 


equations, the electron fluctuation density and particle electrostatic potential are 


én (kw) = a [o(kw) + o(kw)/elx (kw) 
5 ape ¥ 5 Eq 25 
-4n 
k = re.6n 
o(kyw) = pr Zesén, 
The ponderomotive potential, ¥(k,w) which drives the low frequency electron 
fluctuations is 
= e@¢ ; : 
ye mu? |Eo+E- Eoeeee| Eq@2.6 
Substituting Equations 2.5 and 2.6 into 2.4 and eliminating On, in 2.4 finally 
yields the general dispersion relation 
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where 
euler ce ene 
D. c<ke x 2K koe 20, w 
a = 1 + X_(w, ok, ) + Xz (w, ok,) 
nd vi. = (eE /mw)'/2 the elect i locity. The ion VI 
a Os of e%o , the electron quiver velocity. The ion Vlasov 
equation was also used to eliminate én. in 2.4 and ion ponderomotive force terms 


were ignored. The usual approach is to disregard the nonresonant antiStokes components 
(w, sk.) and for Brillouin (and Raman) scattering, the electrostatic wave is 
coupled to two electromagnetic waves. The reduced dispersion relation is then 
1 ] k*v,.?cos*0 
Sea Selita a bad Eq. 2.8 
Trt X D_ 
where © is the angle between the incident and scattered electromagnetic wave 


vectors. From D in Equation 2.7, with |w| << W6 , which is certainly 


true for an ion acoustic wave, Equation 2.8 can be conveniently rewritten for Brillouin 


sackscatter ( 0= TT ). 
2 iz = 
w = Aw + “pe Vos (1 € x;) Aw = (k-k) ¥ k2/2) c*/w, Eq 29 


When the real part of the second term on the RHS of Equation 2.9 is much less 
than W , the real frequency component is = Aw since k = 2k, 
which can be shown from the weak coupling dispersion relations. The electrostatic wave 


phase velocity is 
Re(w)/Re(k) = (c. + u) 


where U_ is the plasma flow speed and Cs in the weakly coupled instability is 
given by Equation 2.2. However, in general, C. may be strongly dependent upon 
the electromagnetic field amplitudes if the ponderomotive force is comparable to the 


free ion acoustic oscillation restoring force (ie. wv/e = 0(¢) in Equation 2.5). In 
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the limit of the fluid description, Mee> w/k o> \V, , the particle 

Lae lh Cree rs Oe a 
susceptibilities are x; ~ “Yj /w Dacha (kao = I/Ap) _ Upon 
substituting into the Brillouin dispersion relation for strong coupling where 
KV Wp; /Bug >> "or , the solution for temporal growth is easily 


shown to be 
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Eq. 2.10 


This.is often called the reactive quasimode solution [16] The spatial problem of stronger 
wave damping will be dealt with later. 

Temporal growth with strong dissipation of the ion acoustic wave through linear 
ion Landau damping occurs for (27/05) = 0(1) and the strongly coupled 
fluid dispersion relation is no longer valid. The ion acoustic wave is shifted to lower 
frequency and kinetic susceptibilities must be used in Equation 2.9. An iterative calculation 
method was used to evaluate complex frequencies for a given real wavenumber. These 
results describe the initial growth of the instability. Equation 2.9 is normalized and 
separated into its real and imaginary parts which describe the real frequency and growth 


respectively of the instability 
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(2k )*K? Eq. 2.11 
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where K = k/2k, a yee w/2koV, . The ion susceptibility is given by 
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where 
Z = 2/K 


W(z) = z2r(z/v2)//2 + 1 


with T the usual plasma dispersion function (Fried and Conte [17]). Numerical 


computation of Tez) was accomplished with a continued fraction method in 
the region ly| ort where 2 tty and otherwise, a convergent 
series for W(z) (Ichimaru [18]) was used for ml eey cae | . These 


calculation routines are summarized in Appendix A. Double precision calculations were 
necessary for W(z) as it was calculated with a large term alternating series. Initial 
values of Q and S&. were calculated from Equation 2.10 and substituted into 
the RHS of Equation 2.11. The 2, and Q. calculated from Equation 2.11 were 
then substituted back into the RHS. This iterative procedure with the appropriately chosen 


initial values rapidly converged to the solution of the dispersion relation. 


Figure 2.1 shows the solution of the dispersion relation under conditions of 


moderate Landau damping ( ZT O/T; = 6 ) and strong Landau damping ( 
ZT A/T; = 1.5 ). For each graph, the strongly coupled limit 
05/2, e/3 and the weak coupling QaVZT O/T. +3 are indicated by solid and 


dashed lines respectively. It is apparent from these plots that strong coupling dominates 
in the parameter space of interest as the 2./2), = V3 line traces the peak 
growth almost exactly in both cases. As the pump intensity increases, both the instability 
growth rate and its bandwidth increase. It must be remembered that these are exponential 
growth rates and the experimentally observed relative amplitudes of the different 
frequency waves are modified by this growth. It is also clear that for decreasing pump 
intensity or increasing ion Landau damping, the real frequency of the ion wave 
approaches the weakly coupled frequency. This is expected because as the driving 
ponderomotive force decreases, the quasimode approaches the normal eigenmode of 
the free plasma. Also, with increased ion Landau damping the instabilitiy growth is 


diminished and the interaction waves are weakly coupled. 
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Figure 2.1 Temporal Growth Rates of Stimulated Brillouin Scattering 
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Even with strong temporal coupling, the real wavenumber does not deviate much 


from 2k, . Figure 2.1 was calculated for k = 2k, endow ai! 


From Equation 2.9 


dw/w. = v1 - n/n. (c/(c. + u))(1 - K) Eg. 2.13 
This equation shows that Aw/w. is = c/C. times that of | Ak/k 


and that the real frequency is strongly affected by the real wavenumber of the ion 
acoustic wave. In a weak electromagnetic pump field Aw = We and from the 
dispersion relations for the electrostatic wave and electromagnetic waves and also from 


wavenumber and frequency matching 


k = 2k - 2c /cvl - n/n.) i. Sal) 


We = ©. on Jee 
eC “o.0mne 


Eq. 2.14 


From Equation 2.14 one sees that even for strong coupling with significant departures 
from the weak coupling ion sound speed, the wavenumber for significant growth is 
always close to 2k 

After several growth periods 0(1/Im(w)) , nonlinear departures from 
this temporal growth and convection saturate the wave amplitudes. Different possible 
nonlinear effects will be explored more fully in the next two sections, but in any case, 
saturation can occur and temporal growth cease although spatial gain may still be 
sufficient to permit convective growth of the instability. From Figure 2.1, again, one 
notices that in the limit of strong damping on the ion waves, the real frequency moves to 
the weak coupling result as the imaginary frequency diminishes. The coupled mode 


equations given in the next section will be more convenient for describing spatial growth 


when Im(k)/Re(k) << 1 and  Im(w)/Re(w) << 1 . However, notice 
from Equation 2.11 that a, = 0 when Im(x, ) = 0 which provides a 
criterion for spatial wave growth A zero of Im(x. ) is approximately 

Fat) tS Im(k) = Re(k) = 2k, and it is easy to 
show w/Re(k) = 2v. . This regime corresponds to Compton scattering 


where the pump wave is scattered from the heavily damped ion wave [19]. Thus after a 
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few temporal growth periods the Brillouin scattering can evolve to strong convective 
growth. A more exact analysis of Equation 2.8 by Forslund et al [10] shows that this 
strong spatial coupling occurs for WON > Hck , a normally 
difficult regime to achieve experimentally, and the broad wavenumber spectrum is 


centered around 2k, 


2.3 Coupled Mode Equations 


2.3.1 General System of Equations 

The results of the last section showed the general dispersion relation and nominal 
behaviour for parametric instability in a homogeneous plasma. The spatial distribution and 
amplitudes of the various waves were not derived since the dispersion relation only 
provides complex %  , k without consideration for the nonlinear effects required 
to saturate wave growth. From Maxwell's equations and the 2 component plasma fluid 
equations, a set of coupled second order partial differential equations can be derived 
which describes the spatial and temporal evolution of the pump field, scattered 
electromagnetic wave and ion acoustic fluctuations [20]. Apart from being restricted to a 
linear fluid description of the Brillouin instability, the equations are quite general. 
However, considerable simplification is possible if WKB approximations to the mode 


amplitudes are invoked, Le. 


E(xst) =e (x,t) exp(i(k sr + wt + Sok +dr)) CeCe 
Bat x tee Caltex 0 (CK cia cue mace ete. 
av . 
n(x t) = ng(x,t) exp(i(k er + wt + sok sdr)) + c.c. Eg 2.15 
c.c. = complex conjugate 
where Gentes e and ne are assumed to be slowly varying quantities in time and 


space as compared to the corresponding frequency and wavenumber of the mode. The 
integral expression is for slowly varying phase changes due to wave dispersion. 
Substitution of these expressions into the second order differential equations reduces 


them to first order equations describing the field amplitudes of the waves. The three 
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simultaneous coupled equations obtained are 


a) 
2a/at + C3a,/ ox = “®p@ Me /WyNgra_exp(-iSOk,-Ak_-Ak, dx)/4 


“9% 
9a_/at + c_9a_/ax Woe Ne/woN,-aexP( isdk,-Ak_-Ak. dx)/4 
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where ‘e is the ion wave damping coefficient. In Equation 2.16, the electromagnetic 


field amplitudes have been normalized to 


er: (k_/k)/* eave See aca) En 


where EG is the peak incident amplitude of the pump electromagnetic field. 
Electromagnetic wave damping has not been incorporated in Equation 2.16 since in most 
instances of interest here, inverse bremsstrahlung absorption is weak in the Brillouin 
interaction zone. Phase mismatch between the waves is accounted for through the 
exponential factor on the right hand side of Equation 2.16. Mismatch sources include 
density gradients affecting the electromagnetic waves whereas velocity and temperature 
gradients cause changes in the ion mode frequency. The dispersion relations for these 


waves Clearly show their dependence on these plasma inhomogeneities, i.e. 


k=- - n/n Electromagnetic Waves 
EG. \2aV7. 


= + ° Ion Acoustic Waves 
anne S45 baal aE 


Further reduction of Equation 2.16 can be effected through non-dimensionalizing time 
and space coordinates and renormalizing the density fluctuation. Different time and length 
scaies are appropriate in the two cases of temporal or spatial problem. The complex 
wave amplitudes can then be separated to rewrite the three complex equations into a 


system of six real equations. The resulting system of equations is then 
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(le. |/c)(aa,/at) + (c,/c)(aa,/az) 


({c.|/c)(aa_/at) + (c_/c)(aa_/az) = aa cos(H(z)+¢) 


Ba, + da./at + dac/ar = aa_cos(H(z)+¢) 


-a_a.cos(H(z)+¢) Eq. 2.18 


(le.|/c)(20,/at) + (e,/c)(80,/82) = (a_a./a,)sin(H(c) +o) 
(le,|/c)(ae_/ar) + (c_/c)(30_/8z) = (aja,/a_)sin(H(z) +6) 


a0 fot + 30,/ae = (a.a_/a.)sin(H(c)+4) 


oh ania acy As 
The complex wave amplitudes have been expressed as e on exp(i0.) 
where es Og Gr cS and a, 95 are real quantities. The ion 
wave amplitude renormalization is 

%~% 
os VZk 72k. (9/9) (Ve/Vog) (L/No) 
and the new space and time coordinates are Gs x(y/ vec.) and 

(ae t(y,/yYc/c,) ; 15 is the homogeneous instability growth rate. 

The gradient phase mismatch integral is H(z) = fhk, ial Sooners Ak. dz and 


B is the normalized ion sound wave dissipation term 
B=VTk.77k5) (w/a) (V/Vog) (ye /ug) Eq. 2.19 


Various simplifications to Equation 2.18 can be made for different special cases of 
instability growth. The usual approximation is to assume phase locking for maximum wave 
growth. Here, aah and the cosine factor on the right hand side of the wave 
amplitude equations is maximum so that wave coupling in Brillouin scattering is largest. 
Plasma gradients can still affect phase mismatch through H(z) and limit the 
spatial growth region. Only cases of $ = 0 and simplifications to either temporal 


or spatial problems will be examined here. 
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2.3.2 Spatial Coupled Mode Solutions 

A full solution of the remaining three equations in Equation 2.18 can be obtained 
directly and the time asymptotic spatial distribution of the wave amplitudes can be 
computed. This steady state Brillouin interaction is the solution to the coupled mode 
equations where time derivatives are ignored. The simplified steady state, homogeneous 


plasma equations are 


304/36 =-aa 


da_/9e - a0 (c_/c = -1) 


‘ Eqi2.20 
30./3¢ + Bal = OL 0% 


with boundary conditions 


The first condition states that the pump electromagnetic wave enters the interaction 
region with full intensity at c= 0 while the second is tne reflected 
electromagnetic wave amplitude as it exits the unstable region. The power reflectivity, 
R is normalized to the incident pump intensity. For subsonic plasma flow (see 
Manheimer et al [21] for supersonic flow), the ion acoustic wave experiences spatial 
growth as it convects in the +z direction with initial amplitude ¢' = 0 for 
ae aA) 
Equations 2.20 are numerically integrated and the solutions are plotted in Figure 
2.2 for several boundary conditions and damping parameter 8 to show how the 
spatial wave distribution and amplitude are affected by these parameters. The amplitudes 
OF and a? are plotted since it is the wave intensity which is 
normally detected in an experiment for both the electromagnetic waves (with suitable 
light detectors) and ion acoustic waves (through Thomson scattering diagnostics). The 


spatial invariant a? -a2=1]-R can be identified from the first two 
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Figure 2.2 Solutions to the Spatial Coupled Mode Equations 
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equations of 2.20. There are several interesting features to these plots. First, for the 
same electromagnetic wave reflectivity, the peak ion fluctuation level decreases 
significantly for increased wave damping. The compensating effect is for the unstable 


region to enlarge so that there is sufficient growth of the reflected electromagnetic 


wave to attain the specified R . This is evident upon comparing the 8 = 1 and 
B Fe5 spatial distributions of uit fore k = 0.1 tienmitc 4 
OL vanishes for 8 = | whereas. for 8 = 5 , even at 
eS: > rie is 0.01 which is still too large for most plausible noise 
sources from which O._ may convectively grow in the as . direction. 

A second observation from Figure 2.2 is that for pe Sel , absolute 
instability can exist in steady state since anil can be zero and yet the waves grow. 
Foret R 102 18orand y Bee pt Soe is zero at G = 2.4 and grows to 
0.1 at Se =60 . The ion acoustic wave does not simultaneously vanish but rather 


convects out of the unstable region and is damped. In general, absolute instability can be 
created for Rit<eZ , but with increased ion wave dissipation the instability is 
convective (Forslund et al. [10]). 


Another important observation which will lead to further simplification of the 
2 


coupled equations is related to the ion wave distribution. From & = 0 ae 


grows rapidly to its peak value because both a and O_ are near 
maximum and the coupling on the right hand side of Equation 2.20 is large. As ¢f 
increases, the ion wave amplitude passes through a maximum and decreases with 
declining O . What is significant is that for large wave dissipation %, is 
more slowly varying through the unstable region. This suggests that for 
peSSen , the ion fluctuation can be approximated as constant and the spatial 
derivative in the ion wave coupled equation ignored. With this assumption, an analytic 
expression for the reflectivity from a plasma slab of length L can be derived from 


Equation 2.20 (Tang [22]) 
Eqé2-2)1 


R(1 - R) = e(exp(G(1 - R)) - R) 


where 
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k Lon ni) Le wienord Eq. 2.22 
pe 6, 5 
G=——|ji-— —= —— 
2 
4 n nN sic 2y./w. 
2 2 a . 
and OL ne is the boundary condition representing the scattered wave noise 


source from which the instability convectively grows. The fluctuation level 


Av 
no/Ny is now related to the reflectivity through 


n/N z 0.5(v../v,)2(1/[2y./w.])¥R Eq. 2.23 


The most obvious source of red shifted electromagnetic noise would be reflection from 
a moving plasma boundary. If plasma densities greater than we are created during 
plasma formation, unity reflectivity of the pump electromagnetic wave can occur. Low 
reflectivity from plasma density gradients is also a significant electromagnetic noise 


source. It can be shown (Born and Wolf [23]) that the reflectivity R of a plasma 


boundary with a linear density gradient from n/n. to O over a distance L 
isgivenby R = Take where, 
bier dit ¥ 
r -| — — exp(k, iG N(z) dz) dx 
He Eq. 2.24 
and N(x) = v1 - xn/nb 
For example, R =3.6 x 19074 for n/n. = 0.5 and Kot = 25 


which, as will be seen in the data analysis of Chapter 4, is a reasonable noise level to 
seed convective instability. 

The ion acoustic wave damping Y,_ in Equation 2.22 can be treated as a 
phenomenological damping coefficient. However, if the electron and ion temperatures 
are known for a Maxwellian plasma, the kinetic ion Landau damping coefficient can be 
substituted for the hydrodynamic value resulting in 


wg/2v, = ImE(1 +x,)/(1 + xq + x9] ae 
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This result can be readily shown by comparison of the coupled mode growth rate to the 


temporal growth rate from the kinetic dispersion formula, ie. 


y= oe lal COUPLED MODE 


(woe /8uy (Vo ?/¥,7) Im[(1 + x /0 + X + agit KINETIC 


Offenberger et al. [24] used the random phase approximation (see Tsytovich [25] for a 
thorough explanation of this approach ). to arrive at the same expression for a 
monochromatic pump field. This earlier result demonstrated the equivalence of the 
random phase mode! and the coherent pump model in the limit of strong ion wave 


damping. Substitution of the kinetic expression for the electron and ion susceptibilities 


results in 
We Z1) xX exp-x2 
ee T. [1+(Z2T,/7;) (1-2xexp-x? [exp(t?)dt) J +(2T,/T; )??xexp(-x?) 
fe) 
Eqe2zi2e6 
where 


x = WM{72K,T, (w_ - 04)/|k_ - 


k,| 


Real frequencies and wavenumbers are used to evaluate X since the coupled mode 
equations apply to weak spatial coupling for which Equation 2.26 is being applied to a 
ko ni Ik. and in 
Equation 2.14 was shown to be very close to 2k, . Equations 2.21, 2.22 and 2.26 


steady state solution. The wavenumber difference is | k_ = 


thus give gain as a function of the ion wave frequency shift through Y./W< 
What remains to be done to complete the description of the heavily damped ion 
wave model is to determine the electron and ion temperature ratio and the fraction 
n/n of the ion distribution which is contributing to the wave damping. Two 


_approaches can be taken. The first approach uses a very simple ion heating model 
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proposed by Phillion et al. [26] where the fraction of the pump wave power delivered to 
the ion wave given by the Manley—Rowe relations, W/W, , i§ Converted to 
thermal energy for the entire ion population in the interaction zone. The effective ion 
temperature is calculated through balancing this input power with energy flow out of the 
unstable region as the heated ions stream out at their thermal velocity. This energy 


balance can be written as 


BG wesdd, 


where I is the electromagnetic pump intensity. With We determined from the 


weakly coupled ion acoustic dispersion relation and the above expression, 
ZT O/T; can be determined for a given pump intensity and reflectivity. For 
reflectivity calculations it is more useful to express the pump intensity as the dependent 


variable whereby it is then easy to show 


Vote Wiese 3T.|1/2|-1 
le i cee eaez. eta Eq. 2.28 
Ve no R [ T; dT, 


Substitution of Equation 2.28 into Equation 2.22 together with Equation 2.26 is the ion 


heating model for steady state Brillouin scattering. The result of this substitution is 


Kot n Li/2inef2 ZT, -3/2 Haley? ] eee | 
R(1-R) = eexp|—]1-— —| |— ]+— 
Ne no) Uy 27 J R 2y. N 


Eq 2.29 


A useful extension of this model has been made through introducing the integer 
divisor Nin Equation 2.29. This corresponds to a multiline laser where N__ is the 
number of different equal intensity laser lines comprising the electromagnetic pump. If 
the frequency separation of these laser lines, Aw) is greater than the instability 


bandwidth, Y each line will act as an independent pump [27]. For 


O , 
N> 1 , the reduced pump intensity for each scattering process results in a 
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significant decrease in the total scattered power due to the exponential dependence on 
N 

Relevant solutions for the ion heating model are plotted in Figures 2.3 and 2.4. 
Figure 2.3 shows the pump intensity dependence of the reflectivity at the ion wave 
frequency with maxium gain in Equation 2.29. The curves show the self-consistent 
electron—to-ion temperature ratios due to heating as well as reflectivity calculated for 
various noise levels as a function of pump power. Figure 2.4 is similar to Figure 2.3 
except that a single noise factor eo ies is chosen, but more importantly, the 
reflectivity is shown for two plasma densities and the decreased reflectivity is illustrated 
for multiline pumps | N = 2 and 3). 

Inspection of Figure 2.3 shows that ZT,/T. decreases with the pump 


wave intensity as expected through ion heating. For an oxygen gas target, with effective 


Z = 6, the ions are heated to the electron temperature for 
Z ; 

ev 290 (1) but further increases in the pump intensity are slow to 

decrease ZT /T, since ions are streaming out of the unstable region. In both 


Figures 2.3 and 2.4, an initial decline in the reflectivity is evident for increasing pump 
intensity as ion Landau damping increases. However, at some point the rate at which wave 
dissipation increases is unable to fully counteract the increasing pump and the reflectivity 
R begins to rise again. This transition in the reflectivity function occurs at 
ZT A/T; = 4 . Also the weak influence of the noise level € on the 
reflectivity is evident in Figure 2.3 and the logarithmic dependence can be seen from 
Equation 2.29. From Figure 2.3, a 10° Variation in noise level changes the observed 
reflectivity by approximately a factor of 3. Finally, Figure 2.4 demonstrates the reduced 
scattered light reflectivity which can be obtained by dividing the pump into several equal 
intensity components. For the conditions of this calculation, the reduction factor is close 
to the number of lines participating; for example, the reflectivity decreases by 
approximately a factor of 2 when a 2 line pump is used. The decrease is not exponential 
because a lowered reflectivity moderates the ion heating and consequently ion Landau 
damping decreases. For the modest pump intensities and reflectivities of these 


calculations, the combined effects of pump intensity division and reduced ion wave 


damping allows only a small change in the Brillouin backscatter reflectivity. These results 
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Figure 2.3 SBS Reflectivity with lon Heating Model 


ALEAAES 


2 


‘ 
+ 


| a a i ace ae ae a ae 


: 


ht ttt 


Q 


20 


¢ Ol =? A® OSI =" 


( a /N 90!) ALISNALNI 


ann “ ¢) 
BT 
SANIT 2 21 


a 

oT 

8T 
JNIT | 

Be 


ee ee 
b'OQ =2%uU/su 


ALIAILI91S9a SAS SANTIIL INN 


(€Z) 4 


Figure 2.4 SBS Reflectivity in lon Heating Model with Multiline Pump 
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will be very useful in the data analysis where the laser was tuned to oscillate on more 
than one line. 

An alternative saturation mechanism to wave damping is ion trapping in the 
potential wells of a coherent ion wave. The trapping process has been described 
extensively in the literature (see for example, Ichimaru [18] or Davidson [28]) and will be 
outlined here only in the context of dissipating the ion sound waves driven through 
Brillouin scattering. lons whose thermal velocity closely matches the phase velocity of the 
ion wave can be caught in the wave's electrostatic potential and be accelerated to the 
wave phase velocity. In the ion wave frame, these trapped particles bounce back and 


forth inside a single wavelength well and resonantly exchange energy with the ion sound 


wave. The velocity range, Viv of the trapped particles is related to the ion 
wave potential ¢ = Mave ./2e where 4 is the ion wave potential 
related to the ponderomotive force o« Eo or eles . Thus in the laboratory 
frame, ions whose velocities V are between 

oe = Viv < : sce ae vee Can experience trapping. The bounce 
frequency, Wp , of these trapped particles is Wp = k(eo/M. Wee and 
in times 0(1 /wp) the trapped ions are accelerated to = 


There are two interpretations as to how the trapped ions can moderate the 
Brillouin instability. In one model, the ion wave amplitude can grow until a significant 
number of ions are trapped and therefore become an energy sink to clamp the ion wave 
level. A waterbag ion distribution is used to estimate the maximum n/N, of 
the ion wave. Figure 2.5{a) illustrates the uniform ion velocity distribution of the 
waterbag model and the relative position of C. . From the kinetic dispersion 
calculations of the last section, AG (1 + x. )ek?¢ from which, for 

k*a5 <<. | n/N, 2 ed/KT . Together with the expression for 

the trapping velocity Ver , the trapping range can be expressed as a function of 
n/N, . When CA erae V3KT./M. , ions in the waterbag 
distribution are trapped and there is an abrupt turn on of strong ion wave dissipation. The 


fluctuation amplitude at which this occurs is easily calculated from the above (Dawson et 


al [29]) 
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Figure 2.5 Two Particle Distributions for lon Trapping 
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% 1/2 1/22 
Ne/Ny = 0-5E(1 + 31,/2T,)'/° - (31,/Z7,) [24 Eq 2.30 


This equation, together with the Bragg reflectivity relation 


R = (1/2) (B/N, )°(n/n, )2(L/a)? 
determine the saturated Brillouin reflectivity. For example. if ZT O/T. = 6 
gta and n/n, = 0.5 then fi, /N, = 0.134 from 


Equation 2.30 and the power reflectivity, R is 1.2%. 

An alternative interpretation of ion trapping which also provides a useful 
connection to the ion. heating model can be imagined. In this model, the fraction of ions 
trapped are assumed to be heated to a temperature * a , thereby creating 
a second distribution function on which the ion acoustic wave is linearly Landau damped 
(Manheimer et al. [21]). In any realistic laser-plasma interaction, the potential well 
coherence would not persist through the instability and a distribution of ion energies 
around ZT would result. This would not be unlike ion heating in that the ion isn't 
trapped for many bounce periods but is heated and convects in the direction of the ion 
wave. However, strong heating for only that fraction of particles in the trapping region 
occurs rather than the moderate heating of all ions in the ion heating model. This picture 
compares well with particle simulations done by Estabrook et al. [30] where equivalent 
damping occurred with only a fraction of the ions heated. The ion wave damping is thus 
moderated by that fraction of ions trapped, n,/n . Figure 2.5(b) illustrates 
trapping for a Maxwellian ion distribution. The fraction of trapped ions is easily calculated 


by integrating the Maxwellian over the trapping region and gives 


n,/n = 0.5Lerf(v(ZT,/T.)(1 + Y(26n/n)) - sgn(1 - Y(2én/n)- 


ave 
where exp\ 3 a 
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The effective wave damping is now n,/ney/w. where vane is 
calculated from ion Landau damping with iF = ZT, . The ion wave damping 
coefficient, Y/W. is large ie. O(1) and the steady state, heavy 


damping reflectivity equation is rewritten as 


kK Uew, on on n |1/2{v__2 
BileccR) =ne 6X ose: eee eee 5) (1 -R)} -R 
cee AGHnt Me ve Eq. 2.32 
The normalized ion fluctuation level, n/N, from Equations 2.23 and 2.32 


completes this alternative ion trapping model. This model is not expected to accurately 
describe the instability growth at high reflectivities where the ion fluctuation is large. 
Figure 2.6 shows several solutions to these two equations. For increasing 
electron-to-ion temperature ratios the reflectivity increases at low intensities. This 


results from fewer particles being in the trapping region for increasing phase velocity. 


However, for (vo/ve)* mae , Most ions are trapped for all initial 
ZT A/T; and the damping does not change with further increases in 
(vo</Ve) . From Figure 2.6 one also sees an initial decrease in reflectivity 

as a function of (vio/vQ)? for large ZT /T; . This occurs because the 


rate of increase of trapped particles with increasing pump is sufficient for the damping 


% f 
to outweigh the pump and although nN increases, the reflectivity does not. 


2.3.3 Temporal Coupled Mode Equations 

The steady state parametric instability is relatively straightforward to analyse 
because a system of ordinary differential equations describes the behaviour. When the 
instability’s temporal growth must be dealt with, it is no longer obvious that a similar 
simplification of the coupled mode equations is valid. Obviously, if the wave amplitudes 
are constant throughout an infinite plasma, the spatial derivatives can be ignored and a 
system of three differential equations with only time derivatives describes the 3-wave 
decay. The appropriate temporal, spatial and density fluctuation normalizations to reduce 


Equations 2.16 are 


t= 2v2Z(cko/c Kk.) (V/V .)(1/u, dt = Typ xre ORs ie 
n/N = V(Z7T) (CK 7k.) (p/w 6) Voe/Va)%e 
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Figure 2.6 lon Trapping Model Estimates of SBS Reflectivity 
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Note that the time scale is 0(1/y,) where is is the homogeneous 


plasma growth rate. The temporal equations are thus given by 


i! 

0 
R 
R 


3a / at 


da_/at' 


i] 
R 
R 


Se CE 
3a./8t' + B'a =a.a 


where 


Equation 2.33 is very similar to Equation 2.20 describing the steady state problem 
except in one important respect The sign of the driving term is reversed in the second 


equation and consequently the conservation law from the first two equations of 2.33 is 


ee + ie = ] (initial conditions a, (0) = ] 
a.{0) = 0 . Thus the system electromagnetic energy is constant within the 
approximation W/W, << | 


Weiland and Wilhelmson [31] have examined Equations 2.33 and their conclusions 

will be outlined and expanded upon here. When there is no wave dissipation, 
B' = 0 and the solution to the sytem of equations are Jacobi elliptic functions. 
These functions describe oscillatory exchange of energy between the two 
electromagnetic waves and amplitude modulation of the ion acoustic wave. With ion wave 


dissipation however, the oscillations damp out and the system relaxes to 


os eel and ja _| = 1 ; all the ney has been pansiete to 
the decay electromagnetic wave. Figure 2.7 is a plot of a_ and Ot for 
several B' FOr pie ev0 , the wave amplitudes stabilize to their steady 
state values in times 0(10/y,) with some ringing. Increased dissipation limits 


the ion wave's peak amplitude which after several oscillations, becomes negligible. Once 


ne 0 , there is no parametric coupling of the electromagnetic waves and the 


instability ceases with the end result of the pump wave completely depleted and 
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Figure 2.7 Solutions to the Temporal Coupled Mode Equations 
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3. Experimenta! Methods 


3.1 Outline of the Experimental Methods 

A considerable number of diverse experiments were performed in the course of 
this laser—plasma interaction research. There were two broad classes of experiments; 
firstly those experiments from which laser or plasma parameters were measured and 
secondly, those experiments where the laser-plasma interaction was diagnosed. 
Examples of the first type are laser injection mode locking and interferometric plasma 
electron density measurements while the second class included backscatter light and ruby 
laser Thomson scattering measurements. This range of experiments required a wide 
variety of instrumentation and the purpose of this chapter is to explain this apparatus and 
its implementation. 

The experimental facility consisted of two high power carbon dioxide lasers and 
an evacuated target chamber in which a focussing mirror and gas jet target were located. 
Most of the interaction experiments in laser fusion use either planar or spherical solid 
targets. The interpretation of Brillouin measurements from these experiments is 
complicated by large density and velocity gradients in the underdense coronal plasma. 
Also, solid targets often have limited diagnostic access, particularly for forward 
scattered or transmitted signals. These reasons alone motivate one to use gas targets 
when studying SBS. 

Suitable ports on the target chamber allowed entry of the laser light as well as 
viewing accessibility for the different diagnostic experiments. Electronic equipment for 
synchronizing the firing of the two CO, lasers with the gas target generation was located 
in a Faraday cage screen room. Other control instrumentation and oscilloscopes were 
aiso located in this screen room to reduce electromagnetic interference generated from 
the high voltage laser discharges. When operational, the entire laser—target system could 
be fired and data collected once every 45 seconds. In this time the CO, laser capacitors 
were recharged, the target chamber was evacuated to remove the spent target gas and 
then refilled with a low pressure static helium fill. Of course the data obtained from any 
shot was specific to the experiment being conducted. The different experimental 


methods to be detailed in this chapter include carbon dioxide laser injection mode locking 
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and multiwavelength oscillation, laser focussing and target generation, 
absorber/scintillator measurement of plasma X-ray emission, temporally and spectrally 
resolved infrared light detection, ruby laser interferometry and ruby laser Thomson 


scattering. 


3.2 Lasers 


3.2.1 TEA Carbon Dioxide Laser 

The principal laser used for the laser—plasma interaction studies was a high power 
transverse electric atmospheric (TEA) CO, laser. This TEA laser was designed and 
fabricated at the University of Alberta (Burnett and Offenberger [32]). On firing the laser, 
three single stage Marx voltage doublers supplied 100 kV to the laser discharge 
electrodes and pre-ionizing pins. The active volume pre-ionization was accomplished 
through a series of pin discharges located on either side of the main discharge volume 
and electrically connected in parallel with the Rogowski profiled discharge electrodes. 
The total charge in the pin discharge arcs was limited with 570 pF series capacitors. 
Upon firing the laser Marx banks, arcs between the pre-ionizing pins flooded the active 
volume with photoionizing ultraviolet radiation (Sequin et al.[33] ). The resulting O(10* cm?) 
electron density ensured a uniform electrical discharge between the main discharge 
electrodes, leading to efficient excitation and inversion of the active medium. Six pairs of 
opposing graphite and aluminum Rogowski electrodes ballasted the discharge current to 
further inhibit arcs while delivering electrical energy to the main discharge volume. The 
total gain length of the TEA CO, laser was 3 m and the resonator mirror separation was 4 
m. An unstable resonator design permitted 40 J light energy extraction from the laser in 
a 40 ns FWHM gain switched pulse. Up to 1 GW peak power was delivered in this mode 
of operation. Table 1 lists the principal features of this CO, laser. 

The gas kinetics of the CO,-N,—He mixture in the laser is well known [34] and will 
be briefly explained here. Basically, the nitrogen molecule is vibrationally excited through 
electron impact during the laser current discharge and it couples efficiently to the (00°1) 
state of CO, molecules to resonantly transfer its energy and excite carbon dioxide 


molecular vibrational states. The resultant inverted population of CO, molecules (between 
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Wavelength 


Energy Storage 
Capacitance 


Operating Voltage 
Preionization 

Active Volume 

Active Area 
Electrical Efficiency 
Resonator 


Rear Mirror 
Output Mirror 


Gas Mix (CO,:N,:He) 
Laser Light Energy 


Peak Power 
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10.59 microns 


1.2 microfarad 


55 kv 

50 x 2 side discharge pins 
Ls toa 

3267 Cnix® Joffe 

2.7% 


4m unstable 


12.8 m radius of curvature 
5 em diameter 


-4.96 m radius of curvature 
1 cm diameter 


Ler orD 
HO J 


1 GW , gain switched 
1.7 GW , injection mode locked 


Table 3.1 TEA CO, Laser Design Features 
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00°1 and 10°O) amplifies the infrared light and lasing occurs with an optical resonator 
surrounding the medium. Helium is added to the gas mixture to collisionally de-excite an 
intermediate CO, energy level which would otherwise prevent efficient relaxation of the 
COQ, molecule The helium also increases the thermal transfer to the walls. A helium-—rich 
gas mixture was required to eliminate arcs betweeen the Rogowski profiled electrodes 
and the reduced carbon dioxide concentration did cause a slight decrease in the total 
energy extraction. A numerical solution of the rate equations for level populations and 
photon lifetime for the unstable resonator [35] showed a weak He concentration 
dependence of the energy extraction in the first 40 ns of the gain switched laser output. 
Further details of this model will be explained when injection mode locking is discussed 
in section 3.2.2. 

The temporal modulation in the gain switched laser output was different from the 
simple numerical calculations mentioned above. Strong axial mode beating resulted in an 
irregular intensity modulation of the laser output as illustrated in Figure 3.1. Figures 3. 1(a) 
and 3.1(b) show two oscilloscope traces where the laser power was monitored with 300 
ps resolution and 3.5 ns resolution (for estimating the average laser power). The fast 
power modulations as seen in Fig. 3.1(a) can complicate some of the laser—plasma 
interaction experiments and necessitated constant monitoring of the CO, laser output 
power. Figure 3.1(b) shows a limited “nitrogen tail” in the laser output where =80% of the 
laser light energy was extracted in the 40 ns gain switched pulse. A power calibration of 
the laser detector was straightforward from these oscilloscope traces once the total 
laser energy was known. 

The CO, laser's positive branch unstable resonator optics were designed for low 
output beam divergence and good energy extraction from the active volume. Zeroth 
order calculations from Siegman’s equations [36] were used to determine the radii of 
curvature and mirror sizes for 85% laser output coupling. The calculated values are listed 
in Table 1. The cavity mode structure was not predicted from these calculations but 
requires complex calculations based on repeated diffracting aperture calculations [37]. 
Although these computations were not done here, some useful observations could be 
made. Firstly, the problem with the zero'th order resonator calculation is that the laser 


output is presumed to be a uniform annulus. In fact, the transverse modes of the 
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Figure 3.1 Gain Switched CO, laser Output 
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cylindrical resonator are described by generallized Gauss—Laguerre polynomials [38]. 
Rather than rely upon detailed calculations, these radial intensity modulations in the beam 
profile were approximated as Gaussian weighted sinusoids whose period was 
determined through burn paper patterns. These measurments were needed for the focal 
spot intensity calculations in section 3.3.2. A second important feature of the unstable 
resonator was the diffraction limited laser beam divergence. This led to higher focussed 


laser intensities illuminating the gas target in the laser-plasma experiment. 


3.2.2 Injection Mode-Locking of the TEA CO, Laser 

Two modifications of the 4 m TEA CO, laser were made in the course of the 
experimental work. First, the laser was injection mode-locked to produce a train of short 
duration pulses rather than the usual gain switched lasing. This modification is described 
here. In the second change which is described in section 3.2.4, the laser was forced to 
oscillate at several wavelengths simultaneously in either gain switched or injection 
mode-locked pulses. 

Short duration laser pulses were required for their fast intensity rise in the 
laser—plasma interaction and their high peak intensity as compared to gain switched lasing. 
Short pulses can be obtained through mode-locking a high power laser or amplifying a 
low level light pulse. Preliminary experiments with two pass pulse amplification in the 3 m 
_ amplifier showed severe parasitic oscillations and made it necessary to use some other 
pulse generation technique. This was not too surprising as the high gain module was 
designed for laser oscillation rather than light amplification. Conventional amplifiers 
normally consist of a series of low gain modules which are isolated to prevent parasitic 
generation. 

Injection mode-locking [39] was selected for producing high intensity pulses 
from the 4 m TEA CO, laser. This mode-locking method was preferred for its ease of 
operation and did not require modulating high light fluxes with expensive, easily damaged 
electro-optic components. Figure 3.2 is a schematic of the injection mode-locking 
system constructed in our laboratory. The two laser oscillators were a 1 m unstable 
resonator (0.5 m gain length) TEA CO, laser injection oscillator with 100 MW peak output 


power and the slave oscillator, the 4 m TEA CO, laser previously described. 
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Figure 3.2 Injection Mode Locked CO, laser 
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The injection oscillator was a copy of a single section of the slave oscillator but 
with Brewster angle sodium chloride windows rather than the resonator optics sealing 
the laser box ends. This was done to produce vertically polarized light from the injection 
oscillator, making it suitable for pulse slicing in the front-end electro-optic shutter. A 
laser triggered spark gap (LTSG) provided automatic synchronization of the shutter with 
the injection oscillator and controlled the light pulse width. Upon firing the injection 
oscillator, the LTSG generated an =9 kV, 2 ns charge cable voltage pulse to activate the 
Chromium doped Gallium Arsenide Pockels cell {7 mm x 7 mm x 40 mm). For the 2 ns, 
the injection oscillator output was rotated 90 degrees upon traversing the Pockels cell 
and subsequently analysed by a 3500:1 polarization rejection germanium Brewster stack. 
Low level unpolarised light from. the injection oscillator limited the pulse-to—leakage 
contrast to 200:1. 

The switched out laser pulse was directed into a 2 mm diameter hole drilled into 
the rear mirror of the slave oscillator and diffraction coupled to the laser mode volume. 
The timing between the injection and slave oscillators was electronically controlled 
through two delay generators. Timing jitter of 10-50 ns resulted from shot-to-shot 
variation in the two oscillators’ gain build up but was not critical to achieve mode-locking. 
The light pulse switched out from the injection oscillator output was injected into the 
slave cavity at threshold, ie. before spontaneous emission build up could lead to normal 
gain switching. The subsequent cavity photon growth is described as regenerative 
amplification. The injected light pulse grows in the laser cavity and only a fraction of the 
pulse energy is re-introduced into the cavity after reflection off the output coupler. The 
resultant laser output is a series of short pulses separated by the slave cavity round trip 
transit time. The pulse width is approximately that of the injected pulse and the pulse train 
envelope is similar to the slave oscillator gain switched output. 

Table 3.2 summarizes the significant IML system characteristics. As seen in Table 
3.2, the IML laser peak power was nearly double that delivered by gain switching. The 
pulse-to-low level lasing was highly variable and was a function of both injection timing 
and the injected photon density. Within the 100 ns injection window, the pulse—to—noise 
contrast was typically better than 20:1. Further improvement in the contrast was 


accomplished by inserting a low pressure (=4 Torr-cm) sulfur hexafluoride saturable 
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Injection Oscillator -500 mJ, 40 ns gain switched output 


-Self synchronized laser triggered 
sparkgap to GaAs electro-optic 
modulator 


Injection Characteristics - 10 W total power 


Slave Oscillator 


- 50 ns injection jitter 

- 100 ns injection window 

- 2 to 3 ns FWHM injection pulse 
= O02 TOs Ue) am gain 

- Peak power to 1.7 GW 

- Horizontal Polarization 


- Pulse contrast better than 20:1 


Table 3.2 IML Laser System Performance 
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absorber cell between the injection and slave oscillators to reduce low level light leakage 
from the germanium analyser. IML laser prepulse levels were difficult to control due to 
its sensitivity to injection characteristics. This meant that the CO, laser power had to be 
constantly monitored during laser—plasma experiments. One virtue of this varying output 
power was that pump intensity scaling of experiments was easily done- a large number 


of laser shots would provide data points over a large span of pump intensity. 


3.2.3 IML Laser Model 

A computer model was written to assist in understanding the injection 
mode-locking and to conveniently study the effects of different operating parameters. 
The CO,-N,—-He rate equations used in these computations are similar to Manes et al.[40]. 
A cell model was used whereby the active medium was partitioned into cells of length 

ola , where c is the speed of light and to the injected pulse width. Competition 
with spontaneous emission was ignored and the slave oscillator resonator mode 
structure was assumed to be spatially uniform. Consequently, the model did not predict 
the observed injection window (see for example [41] where a different model is used to 
obtain this) nor the complex spatial variations in gain saturation throughout the active 
medium. The injection pulse was introduced R1 round trips after the start of the slave 
oscillator discharge current and advanced to successive cells each time step. The rate 
equations for each cell were solved with zero photon density in all cells except where 
the pulse currently resided. A fraction of the photon density, 1-R, coupled out of the 
laser cavity each time the pulse reached the output coupling mirror of reflectivity R. The 
remaining photons passed through the cavity twice before a second laser pulse was 
produced. Zero gain cells could also be introduced at either end of the active medium to 
model the effects of changing the cavity length. 

A simplified energy level scheme for the CO,-N, kinetics is given by Andrews et 
al.[35] and was used here. As seen in Figure 3.3, the nitrogen is excited through electron 
impact to level N which resonantly transfers to carbon dioxide level n,, at rate k. Within 
the CO, energy levels, there is stimulated emission of the laser line from n, to n, with 
subsequent decay to ground n, through intermediate energy level n;. Electron impact also 


affects the CO, n, population to some degree. The various rate coefficients used in this 


| - ; ae een 
gS - : | ‘ the et | rae bis oi Te » | 14g = a : 


ad..of Went Te wog jWes! 03 arity fet fia ie 2b 
lea prtiviaw arte to’ eusiy afio sfrarninagxs “ ve 


Yedinun. eD"e/ & —Ssrob vik 


Toi bint wy onfereictet fii Feds OF PSL TM deli” ‘eam, 


there, 5 a 


ry 


2xslOme 142 ° iB tW;. Tans seh © te 4 . etqstig e 1 RE pe yhawe mT "BYRDS. . i C¥A ye i. 6 : 
| ee ahr 
WoMs TS SSrsiv) Di Rive O16 2ifRIGK ier) mol 3 on, ni pady ‘aceiniipe ate ; 


ricws 1 Tg) OT Gerrit 2AM muna ayit i . at weet baa saw Inve ‘ 
Yo hw Baws os1Oe(hh.e | bem a‘ a ae wig ‘cae | phe ¥ io ae nV : 


On. 2OlAno2s1" Wpenines ovata any ‘bos ‘be xing! aa 12 de 
o913 10" bin Iekgarm sk “WN8 gaa mae ieee ornare si 
pasty ai Isbor iis iis a 18 Wy bet sighres ler-goal w Sharir \ 

t tudes iit nomena eg, i straaenay lease si 

evel sat to Mnete att Site aay baud | anyon aha tial 


a 


} i 
Bah d ae an 
= ph ie 


6/67 on gets Srl’ Toss gies eviecoogue, wav 
: ; a ‘ @ % 


Biectw }op3Ns aon Ne Al vi ieigk rioteriy, ayes, vieaah - 
ac} io tus belaueS Hef viene Q ona’ on VS ede 4 * 
| ytivitosiad #¢ Aone pelighion fudituo, are vatisot stg] . 


=) i eeaaae 8 ie 


ZEW a2h0 ise8t Dooase, & Stoted a2 Seat oar “uit 


of myiban! evitse ert +4 bre vefthe.te psoubaieh i 8d Sebbibiu0s atid ome 

at | “tinal ita Saif ‘grigheris toa 
1@ 2werhA yd novig 2 eaiteaid 7 M9, 0 107 amadoe lave! vera: oii 
nowasle Houe wl! Oafioxs ei Ndpamin st es out fa nage 2 oa anert bees ‘ a 
Ain We 9163 7a ,N lgvai abixaib nodisa, or ergtenaat gitar rir | 
AAW cA Ob af Mot, anil seer ort TS rowan’ berehurmnt:. ai. srgitt Alaver 
gals Jonarti NGYseIS |r Iavel vO: sn Sao Gdn yy ono Sia 


ait ni Dezu ainegi tos 918% i ones i 00D a 


Patt - 
a u hr 6 


ée, 


Figure 3.3 Simplified CO,-N, Energy Level Diagram 
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model were calculated for the 73:5:22 He:N,:CO, laser gas mixture. The rate and 


conservation equations used for each cell are: 


dn, /dt = 6yn(n, - fny) +K(Nn, ~ nino) - K13") - ocq(n, - no) 
dN/dt = bon (N, - fN) - K(Nn, - nN) 

dno/dt = ocq(n, - no) ~ Kon, 

dn3/dt = Kono - K3n3 + Ki 3") 


dq/dt = ocq(n, - No) 


n =N Saya = in 
fe) C02 ] 2 3 Eq 31 
Ny i Nye oan 
ec No (4t/7) 01 - t/t) AER a 
= 0 otherwise 
where Nn is the discharge electron density and q the photon density in a cell. 


e 


Nco, is the total CO, molecule density and similarly Nyy, is the total N, density. The rate 


coefficients and cross sections for this laser are: 


= 6 x 107? cm? s7! 

65 ee X 10°78 cm? s~! 

CER ie cm> s7| 

Ky3 = (85(He) + 110(N,) + 365(C0,))760 = 1.12 x 10° om 

Ky = 1.06 x 10°(0.056(He) + 0.46(Np) + (CO,))= 3.07 x geet 
6 ol 


Kz = 760(4000(He) + 40(N.) a 200(C0,)) 202 ¥25" XPT0 


Got xX 107!9 em7* 


f= 1.157 


(He), (No). (CO,) are the fractional helium, nitrogen and 
carbon dioxide gas constituents 


The laser electrical discharge has been described by Burnett and Offenberger [32] who 
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measured a 400 ns current pulse. Electron density calculations [42] from their figures 
indicate peak electron densities in the discharge of 3 x 1033 cm’. 

Equations 3.2 were solved numerically with a second order Runge-Kutta 
integration and the principal conclusions from the modelling are as follows. There is only 
an extremely weak dependence of the peak IML laser output power on injected photon 
density and timing. This was tested with six , 2 ns gain cells comprising the laser and 
energy from 10-7 to 10-* J injected into the slave oscillator at 1 to 18 round trip times 
after the start of the laser discharge current pulse. Spontaneous emission was not 
incorporated in the model since the injection window is not as large as these results 
indicate. Lowered output coupler reflectivity reduced the cavity feedback and moved the 
peak pulses into higher gain above threshold. The outcome was higher peak powers in 
the IML laser pulses with lower coupler reflectivities. A larger resonator gave the same 
effect, with an 8 m cavity giving =40% higher peak powers than for 4 m resonators with 
the same active medium length. This behaviour has been experimentally observed in other 
laboratories [41]. Another very important effect clearly demonstrated in the model was 
the dependence of the peak pulsed power in the injected pulse width. The energy per 
pulse remained approximately invariant to the pulse width, resulting in decreased peak 
powers for increased pulse width. This effect was very noticeable in the laser operation 
where poor electro-optic gating of the injection oscillator output reduced the peak 
intensities by factors of two and more. This made it essential to pay close attention to 
the alignment of the laser triggered spark gap and electro-optic shutter. 

The model predictions compared well with the initial pulses of the observed IML 
laser output. A model result and two measured laser pulse trains are depicted in Figure 
3.4. After the first few pulses, the mode! does not closely resemble the actual IML laser 
output After gain depletion from the first four or five laser pulses, the oscillator gain 
falls below threshold and the spontaneous emission competes with the IML lasing in the 
nitrogen tail. Consequently, the nitrogen tail was reduced to a weakly modulated pulse 
train as significant low level lasing decreased the gain available to the IML laser pulses. 
Figure 3.4(c) shows how SBS fed back into the slave oscillator caused double pulsing in 


the nitrogen tail of the laser output to further limit the late time peak powers. 
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Figure 3.4 Comparison of IML Laser Model to Experimental Results 
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The utility of the IML laser model came in identifying the parametric response of 
the IML laser to injection timing and density, resonator coupling and length and injected 
pulse width. These results revealed the important features to attend to in operating the 
lasers. Since only the first high power pulse was important in the laser-plasma 
experiments, the late time deviations of the model from the real laser performance were 


ignored. 


3.2.4 Multiwavelength Oscillation of the TEA CO, Lasers 

The CO, laser operation described above has been for single line extraction, that 
is, of all the available carbon dioxide lasing transitions, only the P20, 10.59 micron line 
was produced. This was simply a result of P20 having the highest gain and therefore the 
fastest buildup to gain switching. Upon gain switching, the CO, rotational levels 
re-thermalize in =0.2 ns and thus gain depletion occurs for all the possible transitions as 
energy is extracted by the dominant P20 line. Single wavelength oscillation is not always 
desirable. With extremely short laser pulse generation, the thermalization time can be too 
long for effective energy extraction from the active medium. Secondly, in a CO, 
laser—plasma experiment, a multiwavelength laser can reduce stimulated Brillouin 
scattering as was demonstrated in this research. It was also of scientific interest to 
demonstrate that it was possible to obtain multiline injection mode-locking of a high 
power CO, laser. This section will describe the criteria that need to be satisfied to obtain 
multiwavelength laser oscillation and how it is accomplished. 

Two methods have been used in conventional gain switched TEA CO, lasers to 
obtain multiline lasing. Cavity etalon tuning provides line selection and can preferentially 
detune P20 to allow other transitions to lase [43]. The second and more common method 
is to install intracavity absorption cells which selectively absorb the P20 line and inhibit it 
from lasing. Our choice was absorption cells with sulphur hexafluoride (SF,) absorption 
gas. SF, has strong absorption at 10.6 microns as compared to other CO, laser 
transitions and, as a saturable absorber, it conveniently bleaches at laser gain switching to 
minimize laser absorption losses and maximize extracted laser power. Figure 3.2 
illustrates the 7 cm aperture SF, cells installed in both the injection and slave oscillators. 


Whether operated in a gain switched or IML output, the slave oscillator could be 
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operated multiline. By simultaneously tuning the injection and slave oscillators to lase on 
common lines, multiline IML lasing was obtained. This is believed to be the first 
demonstration of multiwavelength injection mode-locking of a high power TEA CO, laser. 

Multiline oscillation of the slave oscillator occurs through the simultaneous growth 
of several laser lines which have close to the same gain and consequently gain switch at 
the same time. For the simplified case of a pulsed laser with a linear gain increase from 


Zero: to"'Ga cm! (averaged over cavity) in time T , the gain switch time 7 for an 


injected (or spontaneous emission) photon density ‘ cm-? is [44): 


T = 1 In M/al +/[2t In(¢_/¢-)/ca] $7, 7ca] Eq. 3.2 


where M is the resonator magnification, de is the saturation photon density ( bs = 
2 x 10'§ cm?) and L is the cavity length. The first term of this equation is the time to 
reach threshold from the beginning of the laser gain and the second term is the time to 
achieve gain switching from threshold. To induce two laser lines in the oscillator output, 
their gain switched times must be approximately equal to ensure that gain depletion from 
one line does not inhibit the second line. This condition applies equally well to gain 


switched and IML lasers. For two lines j and k to successfully participate in lasing then: 
IT, - T,| < cavity transit time 


Inspection of Equation 3.2 for 7 shows the main controlling factors in obtaining 
multiwavelength oscillation. First, 7 is only logarithmically dependent on the spontaneous 
or injected photon density oe . For example, with the CO, laser parameters in Table 1 
and a gain of 0.025 cm", if 10 W of power is injected at one wavelength, a second 
injected line with 4-30 W would satisfy the gain depletion criteria and still lase. Injected 
energies over six orders of magnitude (10° to 10° J) were tested in the IML laser 
program and the gain switch time did not change over more than =120 ns. This time 


range is greater than what was acceptable for simultaneous gain switching but does 


illustrate the validity of approximation in Equation 3.2 . 
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The second observation is the direct dependence 7 has on © -small changes 
in the laser gain for any laser transition could drastically alter its gain switch time. Since 
a depends on both the absorption cell losses and the laser gain, their adjustment is 
very critcal to obtain multiwavelength lasing. The gain buildup time 7 is only dependent on 
the electrical discharge pumping of the active medium and is the same for all laser 
transitions. Overall, 7 was approximately 600 ns to 1000 ns within the normal laser 
operating conditions. This range was much greater than the permissible =14 ns 
difference in gain switching for two laser wavelengths. However, control of a 
through the SF, absorption cell enables 7 for several laser lines to be roughly equal. 

Obtaining multiwavelength oscillation from the gain switched CO, laser was 
straightforward. The intracavity absorption cell was filled with 2-250 Torr-cm of SF, 
and helium buffered to atmospheric pressure. As a result, upon lasing, the 10.59 micron 
P20 line was eliminated and lasing on 2 to 5 lines in the 10.3 micron or 9 micron bands 
took place. The output energy was reduced from the single line output by approximately 
50% to 17 J because of cell absorption losses and reflection losses from the sodium 
chloride windows. An Optical Engineering CO, laser spectrum analyser was used to 
detect the various lines in the laser output. The laser energy distribution in these lines was 
determined with an infrared multichannel analyser located behind a 0.5 m infrared grating 
spectrograph (see section 3.4). 

Injection mode-locking on several lines was achieved through tuning the injection 
oscillator with its absorption cell and those lines common to the two oscillators were in 
the IML laser pulse train. Table 3.3 gives those combinations of equa! power lines which 
were obtained for various amount of intracavity SF, . The wavelength content was highly 
sensitive to the laser gain and SF, pressure as expected from Equation 3.2. Drift was 
observed in the output line content as either characteristic changed during laser operation 
and thus the spectrum had to be monitored continuously. 

As with the gain switched laser, the IML laser performance was reduced in 
multiline operation. Peak tota! multiline powers of 390 MW were observed with 
approximately equal power in the participating lines (up to three lines) The power 
contrast of the laser pulses to the low level lasing between pulses was =9:1. This was 


significantly poorer than the 20:1 contrast in the single line, 10.59 micron lasing because 
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INTRACAVITY SF LINES 
0 Torr cm P, (20) 
<5 Torr cm R, (26) 


Ry (10), Ry 9614) 
R, (18) » R, (20) 
Ro(14), R, 9(18), R, 9(26) 
Ro(14), Ry (16) » Ry 9 (18) 
Ro(14), Rg(16) 
5-25 elornecn RQ(14) 
Ro(14), Ro(16) 
Rg(14), Rg(16), Ro(18) 


250 Torr cm Py(18), Po(22), Ro(12) 


Table 3.3 Some Line Combinations in a Multiwavelength IML Laser 
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a saturable absorber was not available to uniformly eliminate light leakage from the 
electro-optic switch at the injection wavelengths and obtain a high quality injection pulse. 
The total extracted energy in the IML pulse train was =15 J with approximately 50% of 
the energy residing in the laser pulses and the balance in low level lasing. Only horizontal 
output polarization was obtained with the multiline IML laser. The birefringent mylar used 
with the single line IML laser to obtain random polarization was a narrow 10.6 micron 
bandpass filter and therefore could not be installed at the injection aperture of the slave 
oscillator for multiwavelength injection. 

Multiline operation of both the gain switched and injection mode-locked 
oscillators was valuable in experimentally confirming the pump dependence in stimulated 


Brillouin scattering. The details of this experiment will be given in Chapter 4. 


3.3 Target Chamber 

The principal experimental equipment consisted of the TEA CO, laser system and 
experimental! target chamber into which the laser beam was directed and then focussed 
onto the gas jet target. The target chamber was a 70 cm diameter by 35 cm high 
aluminum cylinder with eight optical/general purpose ports around its perimeter. A 5 cm 
thick plexiglass plate on top of the chamber was removeable to permit access to the gas 
target, focussing mirror and diagnostic instrumentation. Electrical connections, gas and 
vacuum lines were provided through suitable ports and feed-throughs on the bottom 
plate of the chamber cylinder. The gas jet target was centrally located in the target 
chamber and a gold coated 10 cm focal length parabolic concave mirror was located to 
intercept the incoming CO, laser beam and focus it onto the gas jet Details of the target 
jet and focussing mirror are given later. 

The interior layout of the target chamber and the position of its ports are shown 
schematically in Figure 3.5. The main optical access windows used were ports 2, 3 and 7. 
The CO, laser light was directed through the 7.5 cm diameter sodium chloride etalon at 
port 2 which also provided a reflected light signal for CO, laser power, energy and 
spectral content monitoring. Stimulated Brillouin backscattering also exited from the 
target chamber through this port. Ports 3 and 7 were used extensively for Ruby laser 


Thomson scattering and interferometry. 
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Figure 3.5 Top View of the Target Chamber 
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In operation, the target chamber was evacuated to =0.01 Torr with a dual stage 
mechanical vacuum pump and then backfilled with 10 Torr of helium (except in special 
cases) for stabilizing the gas flow in the target jet. The gas target reservoir was 
simultaneously re—pressurized to 2000 Torr of oxygen at the time of helium backfilling. 
This target chamber preparation sequence was automated through a specially designed 
microprossessor controller which was activated immediately after CO, laser firing. The 
evacuation and gas filling solenoid valves were operated by the controller and solid state 
pressure transducers for monitoring the gas pressures were interfaced to it through 
standard analog-to-digital converters. The controller software was designed to allow 
control system recovery from noise interference generated by the TEA CO, lasers as 
well as electrical discharges from other experiments. After the gas jet target was 
produced and the TEA CO, lasers fired, the target chamber preparation cycle was 
repeated. A 45 s experiment cycle time was possible with this arrangement although 


longer times were often required for data recording or experimental modifications. 


3.3.1 Oxygen Gas Jet Target 

The focussing mirror and gas jet assembly were located inside the target 
chamber. The gas jet assembly provided a supersonic laminar oxygen gas jet which was 
irradiated with high intensity CO, laser radiation. The laser ionized gas provided the 
medium in which the laser—plasma interaction took place. The construction of the gas jet 
is illustrated in Figure 3.6. In operation, the refill solenoid opened to allow filling the gas 
reservoir to 2000 Torr oxygen and then closed. Upon triggering the gas target , the dual 
solenoid valve was opened for 500 ms and the laminar jet was formed. A piezoelectric 
pressure transducer detected the gas pressure pulse below the jet orifice and gave a 
timing signal for instrumentation in the screen room and for synchronizing the laser firing 
with the gas jet generation. 

The expansion orifice of the target jet was designed as a supersonic nozzle. Gas 
flow became sonic at the throat and subsequently supersonic upon expansion into the 10 
Torr Helium environment. Gas parameters at the nozzle were 58 Torr pressure, 5 x 10 
cm-? particle density and 500 m s* flow speed. Shadowgraphy was used to observe the 


gas jet boundaries for different helium gas fill pressures to compare profiles with its 
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Figure 3.6 Gas Target Design 
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calculated behaviour. The gas flow exhausted into low ambient helium pressure. 

Under these conditions, an underexpanded jet formed over the supersonic nozzle 
and the potential equation for compressible flow was solved through the method of 
characteristics to obtain flow parameters in the jet (Shapiro [45]). lsentropic calculations 
were used and the region of validity extended above the jet only until frictional forces 
between the ambient gas and jet became significant. Figure 3.6(b) will assist in explaining 
the method used in the free jet calculations. First, in region 1, the gas is just leaving the 


nozzle, with a Mach number calculated from the isentropic flow relation: 


] 2 eas cael ie y *'1 


Myer | Z 2(y - 1) 
where At 

is the nozzle throat dimension (flow is sonic at the throat) and M is the Mach number at 
any position of dimension A downstream from the throat. The ratio of specific heats is 

Y which is 1.4 for diatomic molecules. Once the Mach number at the nozzle exit 
is calculated, all other gas parameters such as pressure, density and temperature are 
calculated through the various isentropic flow equations. The gas outside the nozzle is 
freely expanding and the jet boundary conditions are determined from the static gas fill. 
In Figure 3.6(b), the gas flow boundaries are shown with the streamline directions 
changing each time a rarefaction or compression wave is crossed. Lines a and b are 
rarefaction waves. 

More detailed calculations reveal that these waves ( they are the characteristics of 
the flow potential equations) are actually Prandt! Meyer fans with a continuous change in 
streamline flow direction across them. The Mach numbers in regions 2 and 3 are easily 
found as their pressure must equal recent (no pressure difference across 
a free boundary). Once the Mach number of these two regions is found, then the angle of 
the free boundary can be calculated from the Prandti-Meyer equation: 

w(M) = Atan(/M2 - 1) - VUy +1)7(y - T) Atan’(y - 1)(M? - 1)/(y + 7) 


Eq. 3.4 


and 6. = w(M, i w(M,) where 0. is the streamline direction in region 
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3 with respect to the jet axis. The rarefaction waves are lines across which gas 
conditions change. The angle these waves make to the average flow direction (eg 
(0, + 0,)/2 across a) is the average of 0. for the two regions 


separated by the particular line - 
Py A -] 
a = Atan[(M” - 1)°°] Eq. 3.5 


From these equations it is possible to follow the flow-parameters through the free jet as 
rarefaction and compression waves are crossed. Once the change in flow direction is 
calculated, the Prandti-Meyer equations can be used to calculate the Mach number in the 
new region (region 4). The calculations are simplified by recognizing the jet symmetry and 
replacing the center streamline with a solid boundary. A program listing of the calculation 
method is given in Appendix B. 

Two cases of Py and Pu, are calculated as above and shown in Figure Oaks 
where Do is the O, jet stagnation pressure and Pus is the ambient 
helium pressure. The measured throat dimension was 113 microns and the nozzle exit 
was 460 microns. In Figure 3.7 the density n/n. is indicated inside each flow field 
where, 1 is the electron density if the oxygen in the jet is ionized to Z=6 and 

ne is the critical density for the 10.59 micron CO, laser wavelength. Two 
features are seen in these results. First, there is a high density sheath around a low 
density region and finally a high density core of variable extent, depending upon the 
pressure ratio. Secondly, the CO; laser focus was set up such that the lowest density 
region was not necessarily accessed. The density sampled in the P =2000 Torr, 


O 


Pie =10 Torr calculation would provide 0.78 ne and 500 microns further into 


the plasma it would rarify to 0.29 ne . The important conclusion is that there really is 
an optimum focal position for the CO, laser in the neutral jet and it must be carefully 
selected. 

Shadowgraphic measurements of the neutral gas flow boundary compared well 
with the above calculations. A pulsed ruby laser light source was used to illuminate the 


neutral gas target for these results. The dashed line in Figure 3.7 efor: thes 10eaore 


background helium is the observed oxygen jet boundary. The comparison of the 
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PHe= 2 Torr 
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Figure 3.7 n/n. in the Oxygen Jet Flow 
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theoretical and experimental results is excellent although it must be remembered that 
shadowgraphy does not measure the absolute gas density. Interferometry measurements 
of the oxygen density were not made because of the limited extent of the low refractive 
index neutral gas. At least the theoretical calculation is a useful guide to the expected 


oxygen densities under different operating conditions. 


3.3.2 CO, Laser Focussing 

The TEA CO, laser and gas target performance have been treated in the previous 
sections. The remaining principal component to be described is the parabolic mirror used 
to focus the incoming CO, laser light onto the gas jet. The mirror consisted on an 
aluminum substrate coated with copper and gold overcoated for high reflectivity, high 
damage threshold and tarnish resistance. Mirror fabrication was done in the department's 
optics shop as described by Arnold and Offenberger [46]. The mirror was mounted on a 
moveable, tiltable carriage so that the mirror focus could be adjusted to anywhere within 
the neutral gas jet. 

The off-axis parabolic mirror focussing of the CO, laser light was modelled to 
assist in determining the focussed laser intensities. Because of mirror aberrations, CO, 
laser beam divergence and laser beam mode structure, the intensity distribution was 
difficult to measure directly. Also, measurements from a Helium—Neon alignment laser 
could only give a limited estimate of the mirror focussing characteristics because its 
divergence and mode structure were very different from those of the actual CO, laser 
beam. 

The calculation approach is shown in Figure 3.8. Meridional rays are incident upon 
a reflecting surface y=f(x) whose radius of curvature is r = R(x) . The radius of 
curvature could be calculated from the shape of the parabola but since it is calculated 
with a second derivative of f(x), it is extremely sensitive to errors in f. Thus the radius of 
curvature of the parabola was directly measured in a Ronchii test. The incoming CO, beam 
is incident at an angle 6 to the parabola's axis and is intercepted by the parabola's 
reflecting surface. Ray tracing was done by approximating the reflected light beam as a 
bundle of rays. The center ray oe was traced and the distance F along its reflection 


was used to define the origin of a new coordinate system (ey) whose 
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Figure 3.8 Meridional Ray Tracing of the CO, Laser Focus 
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A 


X' axis was parallel to the reflected ray. All other rays of the CO, laser beam were 


a 


traced and their y' intercepts calculated. By weighting the rays with the CO, 
unstable resonator intensity distribution, the intensity profile along the y' axis was 
easily calculated for any position along x! - A further refinement included laser beam 
divergence. Divergence was introduced in the modelling by dividing each ray into several 
equal angle spaced subrays which were then traced. These equal intensity subrays 


spanned the range of incident angles resulting from the laser beam divergence. 


The measured functions f(x) R(x) and =-—« W(x) where 
W(x) was the laser intensity weighting function were: 
6 ’ 
f(x) = 5 a,x’ a, = -0.023095 a, = -0.00959025 
a a, = -0.0939695 a, = -0.001734167 
a. = -0.024770834 ag = -0.00011277778 


R(x) = 19.4 + 0.07547 x Or<x 2.65 

= 19.0421 + 0.2105) x "22654 x < 3.6 
19.4727 a0. 0909) xmas < x =< 5.8 
sin? (4nx/3.7) 


W(x) 


The 0.3 mrad CO, laser beam divergence was calculated as the diffraction limit of the 
unstable resonator’s annular aperture [23]. Further, focussing diffraction was 
accommodated by selecting a grid dimension close to the mirror's diffraction limit. The 
diffraction limit in the f/3 focussing was approximately 30 microns and so a 20 micron 
grid size was selected. A finer grid would have generated nonphysical high frequency 
spatial intensity modulation in the focal region. 

Figure 3.9 shows the calculated intensity profile through the focal region of the 
parabolic mirror with the above parameters. Sixty to seventy per cent of the incident 
laser power is within a 100 micron diameter spot in an 800 micron waist. This long waist 
is a result of the mirror aberration and as a consequence, a uniform laser intensity 
iiluminates the plasma region of interest during the laser-plasma interaction (at least until 
severe plasma refraction takes place). At the peak intensity (axial position O microns in 
Figure 3.9), the beam FWHM is 60 microns although at this axial position, a significant 


fraction of the incident power is off axis by 150 microns. 
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Figure 3.9 Off Axis Parabola Focus 
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It is important to consider how the CO, laser focus was aligned over the gas jet 
orifice. An alignment Helium—Neon laser for the unstable resonator optics was focussed 
onto a 200 micron diameter microball suspended above the jet orifice. The mirror focus 
was adjusted for maximum occlusion of the He-Ne laser light from the microball. Similar 
focussing model calculations for the He-Ne laser showed maximum occlusion within 200 
microns of the CO, laser focus which is acceptable positioning accuracy for a long beam 
waist. Even for this case in calculating the CO, focus position and intensity distribution, 
once the oxygen target is ionized, refraction and absorption alter the light intensity 
distribution. However, it is still important to estimate the focussing properties of the 


optical system for laser pump calcuations in the laser—plasma interaction. 


3.4 Infrared Light Detection 
Several kinds of infrared light sensitive instruments were used to monitor the TEA 
CO, laser performance and light scattered in the laser—plasma interaction. Four detectors 
were available for time resolved measurements of the infrared light. First, low sensitivity 
coped germanium photon drag detectors were used as monitors for the incident and 
transmitted light from the plasma. Typical responsitivity and risetime were 0.5 
uV/W into 50 2 and < 1 ns respectively. The observed risetimes 
were slower due to amplifier , cable and oscilloscope risetimes so that =2.4 ns signal 
risetimes were normal. Type 7834 Tektronix storage oscilloscopes with 7A19, 50 
2 vertical amplifiers (400 MHz bandwidth) were the principal oscilloscopes used for 
all measurements . Two other detectors were liquid nitrogen cooled, gold doped 
germanium (Ge:Au) photodetectors which were used in applications requiring high 
infrared light sensitivity (0.2 V/W into 50 & ) and medium bandwidth (1 ns intrinsic 
risetime). Various bandpass, cut-on and cut-off infrared interference filters were used in 
front of these detectors during multiwavelength laser experiments. Finally, high speed 
infrared light monitoring was done with 100 ps risetime pyroelectric detectors (5 
uV/W into 50 2 ) connected directly with 1.5 m of coaxial cable to a type 
7104 Tektronix oscilloscope with a 7A29 vertical amplifier. The risetime with this 


combination was better than 360 ps. 
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Time integrated infrared light measurements were done either with a pyroelectric 
ioulemeter for routine laser energy monitoring or with an infrared multichannel analyser 
(IRMA). The IRMA was a PARC type 1204A pyroelectric vidicon with a 10 mm x 10 mm 
scanned area and 33 micron effective channel width. The IRMA controller had an output 
line to synchronize scanning with the CO, laser firing. Accumulated data were displayed 
on both a numerical readout and on an oscilloscope. Spectrally resolved backscattered 
light measurements were made by positioning the IRMA head behind a 0.5 m infrared 
grating spectrograph with a suitable imaging arrangement. Either a 75 lines/mm , 64 x 64 
mm _, 8 micron blazed grating or a 148 lines/mm , 5 micron blazed grating were installed 
in the spectrograph. The lower resolution grating yielded 3.2 A/ chnl dispersion on 
the IRMA with 35 A resolution and the high resolution grating was approximately 
twice as dispersive with twice the resolution, as expected. The low resolution grating 
was normally used to provide a large scanning range and higher grating efficiency. 

The limited resolution in spectral measurements was improved through 
deconvolution of the instrumental slit function from the data points. To accomplish this, 
the spectral data points were first convolved with the measured spectrograph slit 
function. In the first iteration the changes in the data point values are the correction to 
the original data) Several iterations of this calculation procedure were done by the 
"pseudodeconvolution” method [47]. This approach was superior to direct single pass 
matrix inversion to remove the instrumental response because of its low noise sensitivity 
and avoidance of troublesome computational round off errors. The calculation algorithm 
is given in Appendix A. All spectral data in this thesis, including visible light spectra, has 
had the instrumental function deconvolved. A test of the method was done by 
deconvolving an instrumental function from itself. Ideally, the solution should converge to 
a single point In fact the test resulted in a four times improvement of the spectral 
resolution. This was considered very satisfactory performance and justified its routine 


application to spectral data. 
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3.5 Plasma Electron Temperature Measurement 

The infrared detectors described above provided information on the CO, laser 
light and light backscattered from the laser-plasma interaction. To complete the 
characterization of the interaction, the plasma electron temperature, ip , electron 
density N, and electron density gradient ne were determined . The electron 
temperature measurement is described here. 

Once the oxygen plasma was created through photoionization and cascade 
ionization in the gas jet target, the high temperature plasma emitted bremsstrahlung and 
liné radiation in all directions. The X-ray spectra of this emitted radiation was monitored 
through a very simple thin absorbing foil/scintillator arrangement from which the plasma 
electron temperature was calculated. A detection head consisting of a 5 mm aperture 
drilled in a shielding aluminum housing was located approximately 4 cm away from the 
plasma. Various thicknesses (1.75 - 12 micron) of aluminum foil were placed in this 
aperture, behind which was located an NE102A plastic scintillator. The blue light 
scintillations from the NE102A material were coupled to an RCA 8645 photomultiplier 
through a low loss 6 mm x 30 cm flexible fiber optic bundie. The calculated response 
time of this arrangement including oscilloscope, cabling and amplifiers was 3.7 ns (1500 
V_ on photomultiplier). The relative photomultiplier signal as a function of aluminum foil 
thickness was compared to theoretical calculations for a Maxwellian electron distribution 
in an oxygen plasma [48] to provide electron temperature estimates. In all cases, the 


electron temperature was much less than 1 keV. 


3.6 Plasma Electron Density Measurements 

The plasma electron density was measured through ruby laser interferometry. The 
diagnostic method provided a 400 ps exposure single frame interferogram whose fringe 
deviations were Abel transform inverted to calculate the radial and axial electron density 
distribution. Figure 3.10 shows the interferometer in place for the laser—plasma 
experiment. A pulsed ruby laser light source illuminated the plasma and the Fresnel 
mirrors were slightly angled to each other to produce fringes. This arrangement was 
essentially a Fresnel mirror division of wavefront interferometer with a telescope to 


image the plasma region in front of the recording camera lens. This interferometer was 
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Figure 3.10 Ruby Laser Interferometer 
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selected because the optical paths of the reference and test sections were nearly 
coincident, allowing simple beam handling. Also, maximum fringe brightness was achieved 
as Compared to the 50% intensity reduction with division of amplitude interferometers 
such as the Mach-Zender interferometer. 

The interferometer light source was a Korad K1Q ruby laser which produced 20 
ns FWHM, 200 mJ pulses of 6943 A light. An optical shutter consisting of crossed 
Glan-Thomson polarizers, a KD*¥P Pockels cell and a laser triggered spark gap, sliced a 
400 ps pulse from the ruby laser output. The laser triggered spark gap voltage and 
optical alignment were critical in achieving maximum energy and contrast in the pulse. Full 
wave voltage was applied to the KD*P cell in a fast rising step from the spark gap so 
shutter transmission occurred only on the leading edge as the Pockels voltage passed 
through its half wave voltage. The shutter output was expanded with a convex/concave 
lens telescope and directed into the target chamber so that one half of the beam 
intercepted the plasma region and the other half was an unperturbed reference beam. A 
two-lens telescope then directed the ruby laser light onto the Fresnel mirrors and 
localized the fringes at the plasma. An interference filter in the telescope cut out the 
plasma bremsstrahlung emission. The optical configuration of the interferometer was 
designed to optimize film exposure, image magnification on the film and fringe resolution. 
Typical figures for this setup were 1 uJ/em® energy onto the film, 50x 
magnification and 25-50 micron fringe separation. More information on the 
interferometer design is given in Appendix C. 

The ruby laser was electronically synchronized to the TEA CO, laser firing with 
approximately 20 ns jitter. The main contribution to jitter was in the CO, laser gain switch 
time. Timing of the interferometer exposure to the plasma formation was done through 
simultaneous monitoring of the 400 ps ruby pulse on a high speed Optics Technology 
photodiode and infrared backscattered light on a pyroelectric detector. In most 
experiments the ruby laser was timed to illuminate the plasma during stimulated Brillouin 
scattering from the CO, laser — plasma interaction. 

Interferograms from this diagnostic method were photographically enlarged and 
then digitized on a Hewlett Packard microcomputer. Data smoothing and an Abel 


transform technique were used to unfold the refractive index distribution through the 
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2 2 
plasma. The refractive index formula Mele to Ne /e My, , where 
u is the refractive index and Wr is the ruby laser frequency, was 
used to convert the refractive index to electron density. The computational algorithm for 


the Abel transform is given in Appendix A. 


3.7 Ruby Laser Thomson Scattering 

Aside from infrared light measurements, ruby laser Thomson scattering [49] 
provided the most versatile diagnostic method for probing the plasma instabilities 
produced in the laser-plasma interaction. Both thermal and nonthermal electron and ion 
fluctuations can be measured with this technique through appropriate selection of the 
incident and collected light directions. In the experiments done in this research, only the 
ion fluctuations driven by the Brillouin instability were of direct interest. However, some 
general comments on the Thomson scattering diagnostic method will be made first. 

Thomson scattering takes place when a light beam impinges upon a region with 
free electrons such as in a plasma. The electrons are perturbed by the electromagnetic 
field and re-radiate light due to their acceleration. If these electrons are already in 
motion, the incident light is Doppler shifted in the electron’s frame as is the re-radiated 
Thomson scattered light in the laboratory frame. The outcome is a double Doppler 
frequency shifting of the Thomson scattered light. Since the electron scattering may be 
due to either uncorrelated or collective fluctuations, the Thomson scattering method can 
probe both types of plasma fluctuations. The difference wavevector of the incident 
probe beam and Thomson scattered light determines which is detected. Wavevector 


matching of the two electromagnetic waves and the wavevector of the fluctuations 


giving rise to the scattering is given by k. =k, + ke where kis the 
incident wavevector , k. the Thomson scattered wavevector and ke is the 
fluctuation wave vector. If [kern| | , where = An is the plasma 


Debye length, then the fluctuation scale length being sampled is smaller than the Debye 
shielding scale length. As a result, uncorrelated electron density fluctuations such as 
thermal motion are sampled. In the opposite case where [kerp| Soel ; 
long scale length electron density fluctuations are detected. These electrostatic wave 


fluctuations can be either ion or electron plasma modes and are distinguished by the 
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frequency shift of the scattered light. The phase velocity of the particular wave connects 
the frequency shift Aw of the Thomson scattered light and the fluctuation 
wave number, v= dw/ | ke| 

Simple geometrical arguments assist in interpreting Thomson scattering 
experiments. Figure 3.11 illustrates the wavevector geometry for Thomson scattering 
where 0 is the angle between the incident and scattered light. With electromagnetic 


waves whose frequency is much larger than the electron plasma frequency, 


fC [k.| and so [kel a7 k,| sino/2 , or for small 


Ike | 
0 such as for Thomson scattering from Brillouin ion fluctuations, 


[kel = |k te . The associated wavelength shift of the scattered light is 


“(al 
Se 
TT 


(2. v/c)sino/2 where V is the component of the electron 
velocity in the ke direction. For low velocity ion motion, the wavelength shifts 
can be very small, often less than one angstrom. 

The plasma fluctuation levels associated with the sampled wavenumber determine 
the scattered light level. The scattering cross section do written in terms of the 
scattering form factor S(w,k) is do = o7S(w.k)dude _ Here dw 
and d2 are the frequency interval and light collection solid angle respectively and 

oF is the Thomson cross section,8.85 x 10-7 cm. Integrating S(w,k) over 
all a) gives the total Thomson scattering for a selected k_ . Calculations 


for thermal plasma fluctuations [49] show: 


2 


S(k) = Sw ek) dave 1/7 (1 eee Wee CN le nn AA Zot /T.)] 


a = 1/kAp ese) Gus 


where the first term is the electron mode contribution and the second is the ion mode 
contribution. When KAn es , the thermal Thomson scattering is primarily 
due to ion mode fluctuations; conversely, electron thermal motions dominate scattering 
when kan >>] . When nonthermal levels of fluctuations are present from the 
CO, laser-plasma interaction, the form factor can be greatly enhanced and large 


Thomson scattered signals can be detected. It is in this context that Thomson scattering 
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Figure 3.11 Thomson Scattering Wave Vectors in an SBS Experiment 
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was applied to the study of the Brillouin induced ion fluctuations in the CO, laser-plasma 
experiment. 

The Q-switched ruby laser described in the interferometry experiment and a 
scattered light detector were the main instrumentation used in Thomson scattering. A 
simple lens focussed the ruby laser light onto the plasma and a second lens was 
positioned for Thomson scattered light collection. Figure 3.12 illustrates the equipment 
layout, including the beam dump and stray light shield which were necessary to prevent 
swamping the Thomson scattered light signal with stray ruby light. Precise alignment of 
the ruby laser and collection optics was done with the ruby laser alignment He-Ne laser 
illuminating a 200 micron diameter plastic microball positioned at the CO, laser focus. The 
He-Ne laser light specularly scattered off the microball and the light intercepted by the 
Thomson scattering collection lens could be seen by eye, allowing relatively easy 
positioning of the light detection equipment. 

The ruby laser output was monitored by the Optics Technology photodiode ( 

< 0.5 ns risetime) for scattered signal normalization and timing the laser with 
stimulated Brillouin backscattering. The 20 ns FWHM ruby laser pulse made timing with 
the short duration Brillouin event relatively straightforward although fine structure (= 150 
ps period ) was ocassionally detected with high speed instrumentation. With a 7 mm ruby 
cavity aperture, 0.5 J pulses with 0.8 mrad half energy beam divergence were produced. 
Stops inserted in the ruby laser path to eliminate high divergence structure in the laser 
output lowered the amount of energy reaching the plasma to = 160 mJ. 

Three light detectors were used in the Thomson scattering experiments. One 
detector was an RCA 8645 photomultiplier which permitted absolute scattered light level 
measurements through a blackbody calibration procedure. With the modest 3.7 ns 
detection risetime, only time averaged density fluctuation levels were calculated from the 
amount of Thomson scattered light. The second Thomson scattered light detector was a 
PARC type 1205A Optical multichannel analyser (OMA). The OMA target head consisted 
of a silicon intensified target vidicon with an S-20 photocathode. The target plate was a 
linear array of 500, 25 micron wide channels each with =10-20 photons per count 
sensitivity. Thomson scattered ruby light was focussed into a Perkin Elmer E-1 0.58 m 


grating spectrograph whose output was imaged onto the OMA target plate. The 
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Figure 3.12 Thomson Scattering Experiment. Streak Camera Configurations are Shown. 
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.04 A per channel dispersion and the 0.4 A resolution of this spectrograph/OMA 
setup limited its usefulness in the Brillouin experiments. 

The most useful detector for Thomson scattered ruby light measurements was a 
high speed streak camera. The scattering system is illustrated in Figure 3.13. The streak 
camera consisted of a Hammamatsu C973 temporal disperser, C1000-18 high gain 
vidicon camera and a C1330 temporal analyser. Ten picosecond resolution was possible 
with a 30 micron entrance slit and high light gain ( < 6x 103) was achieved with a 
microchannel plate in the streak tube. Streak camera information was digitized by the 
temporal analyser system and displayed on a video monitor. The display field consisted of 
64 windows across the 3 mm streak camera entrance slit and 256 time channels. Various 
data routines of the temporal analyser permitted graphical displays of the streak camera 
counts. Consequently, the digitized streak camera field was three dimensional as time, 
space and intensity co-ordinates were displayed. The space coordinate was also 
calibrated for wavelength or wavenumber of the Thomson scattered light as different 
imaging methods were used. A fiducial was placed to one side of the streak camera 
entrance slit to permit scattered light normalization and timing to the ruby laser output. A 
50 micron diameter core x 30 cm graded index silica fiber (1275 MHz km bandwidth) 
whose free end was illuminated with a fraction of the ruby light was used for this 
purpose. 

As mentioned above, various imaging techniques were used to illuminate the 
streak camera's entrance slit. Spatial imaging was done with a magnifying relay lens 
focussing the collected Thomson scattered light onto the camera slit. This arrangement 
allowed one to follow the temporal evolution of the ion mode spatial distribution. Better 
than 10 micron resolution referred to the plasma source was achieved with the f/10 
focussing onto the entrance slit. In these spatial measurements, it was important to 
ensure that the ruby laser beam uniformly illuminated the whole scattering region. This 
required a slight defocussing of the ruby probe beam at the CO, laser—plasma interaction 
volume. Figure 3.12 illustrates the optical configurations for spatial measurements as well 
as wavenumber and frequency spectra of the scattered light. Defocussed illumination on 
the streak camera slit permitted wavenumber discrimination of the Thomson scattered 


light. The simpler approach with a single lens collimating the scattered light and projecting 
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Figure 3.13 Main Components for Streak Camera Thomson Scattering 
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it onto the camera slit was not implemented as it required a small aperture lens close to 
the CO, laser focus. Both smal! aperture diffraction and obstruction of the CO, laser 
beam by a close proximity lens made this direct solution impractical. Ray matrix 
calculations for the optical components and path in Figure 3.12 show how in this 
wavenumber discrimination method, a poor choice in optical design could introduce a 


source position dependence in the streak camera illumination. In ray matrix notation: 


"out | a : | Min 


e D 


fi out 


where Yr and Yr are the source and detected ray distances from the optic 


in out 
axis and ”” denotes their slope. The ray matrix coefficients A,B,C and D are calculated 


from the product of each ray matrix along the optical path. To eliminate the effects of 


source position on Yr the A coefficient must be zero which imposes the criteria: 


out 


xrerf hlieemi AGer- afi.aa Eq 38 

2 2 1 
. where t is the collection lens focal length, fo is the imaging lens focal length 
and =L is their separation. The required distance from the imaging lens to the camera 
slitis xX . The collection optics configuration was set up to satisfy this condition. 


f = 10cm, f, =30cem ‘LC © = 350 cm X =33 cm). 

For the small Thomson scattering angle in this experiment, the ion fluctuation 
wavenumber was directly proportional to ©  , the scattering angle. Typical 0 
values expected for Thomson scattering from the Brillouin fluctuations were six to eight 
degrees, depending upon the plasma electron density and wavenumber bandwidth. These 
small angles required effective ruby laser beam spatial filtering and dumping to prevent 
stray light interference. After taking these precautions, Thomson scattering was 
positively indicated as ruby light was detected only when a scattering plasma was 
present. Light refraction from the plasma was not important at these angles as was 
estimated through an approximation of the eikonal equation and the known refractive 


index profile measured through ruby laser interferometry. 
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Finally, in exactly the same manner as with the OMA, the streak camera was also 
located behind the E-1 grating spectrograph for scattered light spectral measurements. 
Time broadening across the diffraction grating [50] limited the ultimate temporal 
resolution capabilities. This temporal smearing is the result of path length differences 


across the grating for a given wavelength, and is given by 
At = mNA/c Eq. 3.9 


where At is the time broadening, N is the number of grating lines intercepting the 

r wavelength light and m is the diffraction order. Time broadening, At > 100 
ps is expected from the E-1 spectrograph. Although this temporal resolution was 
sufficient for all but the fastest modulations, the spectral resolution still limited the 
overall usefulness of these measurements. 

All of these streak camera diagnostic configurations were used to completely 
characterize the ruby Thomson scattered light. Subsequent calculations from these 
measurements provided important insight concerning the ion fluctuation content of the 
plasma and, in particular, those fluctuations driven by the Brillouin instabililty during the 


CO, laser-—plasma interaction. 
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4. Results 


4.1 Introductory Remarks on the Experimental Results 

A variety of experiments are presented and analysed in this chapter. The results 
are examined in an attempt to first characterize the laser and plasma conditions under 
which the light scattering occurs and then to explain the nature of the stimulated Brillouin 
instability. The chapter is divided into three sections. In the first section, the plasma 
‘electron density and temperature are discussed, followed by a brief recapitulation of the 
carbon dioxide laser behaviour both when gain switched and injection mode-locked. The 
next section presents preliminary measurements of the Brillouin backscattered light from 
which certain model predictions are validated. Finally, in the last section, more detailed 
investigations into the ion mode behaviour are interpreted. These are based exclusively on 
measurements of ruby laser Thomson light scattering from ion acoustic waves driven by 
the Brillouin interaction. Cumulatively, this chapter identifies the laser—plasma processes 


which determine the space and time evolution of the Brillouin scattering. 


4.2 Interferometric Measurements of the Plasma Density 

Several interferometric measurements of plasma electron density were 
performed with the apparatus described in Chapter 3. Not only was the average electron 
density measured, but the axial and radial profiles were computed through Abel inversion 
of the interferogram. These density measurements represent data which were crucial for 
determining the stimulated Brillouin growth rate and investigating whether plasma 
inhomogeneities or interaction nonlinearities in a homogeneous plasma control the 
Brillouin gain. 

Three distinct laser-target conditions were examined. Two test conditions were 
with a gain switched CO, laser irradiating the oxygen gas target and a third with the 
injection mode-locked laser. The oxygen stagnation pressure : helium chamber pressure 
ratio was either 2000 Torr : 10 Torr or 2000 Torr : 0 Torr for the gain switched carbon 
dioxide laser while a 10 Torr helium background was always maintained for the IML laser 


shots. The two conditions with a 10 Torr helium background were almost exclusively 
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used for further laser-plasma studies as it assured reproducible underdense plasmas 
even though there was a penalty of neutral gas density gradients within the target which 
could restrict the extent of uniform plasma. 

Figure 4.1 shows an interferogram for a plasma within 1 ns of the start of the 
Brillouin scattering event A 700 MW peak power CO, laser irradiated the helium 
Stabilized oxygen target. In the same figure, the electron density is plotted for the first 
few fringes from the leading edge of the plasma The axial density is well behaved in a 
100 micron diameter cylinder within which the average electron density is 4.3 x 10?* cm? 
(excluding the first fringe at the ionization front) and the density gradient scale length, 

[(dn/dx)/n]! =600 microns. The formed plasma extends 1.5 mm axially 
with a boundary approximating that of the calculated laser focus though wider as a result 
of hydrodynamical expansion and photoionization from the intense line and 
bremsstrahlung emission from the laser heated plasma Also, as seen in_ the 
interferogram, the back end of the plasma is shaped like a diverging lens( yw < 1 for 
plasma ) for which beam refraction can expand the laser spot size over its vacuum 
dimensions. A computer simulation with CASTOR [51] supported this idea though 
unfortunately the simulation did not include photoionization. Figure 4.2 is a second 
electron density profile only now measured 3 ns after the stimulated Brillouin scattering 


began. The average axial electron density is =0.7 Nn , which is considerably higher 


G 
than that in Figure 4.1. A combination of Abel inversion error ( OOD Ne ) and 
shot variation could account for the density change. From many interferograms, the mean 
electron density was calculated to be =0.5 ae in the time interval 3 to 8 ns from the 
start of the laser—plasma interaction. 

When the gas jet was operated with O Torr helium, similar electron densities and 
radial profiles were obtained. At first sight, this might seem surprising because of the 
extreme underexpansion in the jet flow. However, the flow Mach number was low and an 
intermediate density core to the gas flow persisted. Plasmas formed under these 
conditions were generally of longer axial extent and the ionization front moved axially 
=38 x 10 § cm/s. With 10 Torr helium background, on the other hand, the ionization 


front moved at a speed < 1x 106 cm/s during the Brillouin instability .This lower 


speed is thought to be due, at least in part, to laser inverse bremsstrahlung absorption in 
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Figure 4.1 Gain Switched Laser Plasma interferogram and Abel Inverted Density Profiles. 
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Figure 4.2 Plasma Electron Density Profiles 3 ns after the start of SBS 
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the low density, low temperature plasma formed away from the interaction zone leading 
to a decrease in laser intensity at the ionization front [52] ( where radiation and 
conduction losses balance the heating resulting in a relatively stationary front). 

The final series of interferograms were made for 2 ns FWHM laser pulses 
irradiating the gas target. Figure 4.3 presents density profiles for the plasma created 
during the Brillouin scattering. Again, the plasma boundary was similar to that expected 
from the laser focussing geometry but there were two striking differences from the 
previous interferograms. First, the ionization was significantly lower than that obtained 
with the gain switched laser. With the short 2 ns laser pulse, the Brillouin scattering was 
prompt because of the higher pump intensities at the time of gas breakdown but 
correspondingly the finite ionization times result in lower ionization. In the gain switched 
laser irradiation case, the Brillouin scattering was delayed =4 ns after gas breakdown. 
This long time explains the increased ionization of the oxygen for the gain switched laser 
experiment. The second feature of the density profiles in Figure 4.3 was the sudden 
density jump from =0.17 net =0.31 Nn. {calculation error + 0:05 

Af in both cases) over 50 microns. This density step was the result of the plasma 
breakdown crossing a Mach line in the neutral gas flow. The electron density gradients 
were not smeared through plasma expansion at this early time of plasma formation (as 
also evidenced through the low electron densities). From this observation it is clear that 
the laser focus position has to be carefully chosen to ensure that the plasma 
characteristics are not complicated by the density structure in the neutral gas flow. 
Positioning the CO, laser focus 300 microns further back towards the laser eliminates 
this problem in the Brillouin interaction zone. Subsequent interferograms consistently 
showed an average electron density no= Os0>n 2 0.05 n as expected for 


Cc Cc 
a low ionization plasma. 
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Figure 4.3 Electron Density Profiles for an IML Laser Irradiated Plasma 


82 


+ i hd i) 
7 : 
he j 
q iy 
f 
1 
i 
a 
‘ 
* 
i 
a 
» 
<= 
- aear sare Sen 
ww 
: { ‘eal _ 
, att Be. 5 e 
D , ” ae mente g 
wen : Seve rt +i 


aes ee . neat — US pa bal A he ale a 
a, an ‘ ee f 
#: Ard a wiki - 
a ets, | 


aril [cpa eni » 


¥ Fi 


‘een 
nee 
é f f mia eich, Y 


= y : , Lair wen 


J 
a? fi 
a alt 
ra \ ¢ 
ay ( am a omer 
dhe 


na rae 2 i A ie 
es al Hi “%, 


Nm, 
ke 4 


83 


4.3 Electron Temperature Measurements 

The second plasma parameter diagnosed was the electron temperature. As 
discussed in the experimental methods, an absorber-scintillator detector arrangement 
was used to obtain relative X-ray signals for varying absorber foil thicknesses. These 
relative signals were compared to those calculated for an oxygen plasma, emitting both 
line and bremsstrahlung radiation at different temperatures and detected through various 
thickness of aluminum foil [48]. 

Electron temperature measurements were done for plasmas created with a gain 
switched laser and the injection mode-locked laser. Figure 4.4 is a plot of the measured 
X-ray signal from the gain switched laser irradiated plasma for several aluminum foil 
thicknesses. For comparison, these experimental results are superimposed on theoretical 
relative X-ray signal curves. The best fitting theoretical curve gives an electron 
temperature estimate of 120 + 20 eV. The accuracy of the temperature estimate was 
limited by the shot-to-shot variation in the X-ray signal at any foil thickness. There was 
no significant temperature change over the range of 300 MW to 1 GW incident laser 
power. With the IML laser, similar X-ray measurements gave Te = 150 + 20 eV for 
500 MW to 1 GW incident laser power. 

All measurements were time integrated by the 4 ns detection risetime and 
consequently the temporal evolution of the X-rays was not examined in detail. However, 
it was observed that the < 20 ns signal pulse duration was shorter than the =40 ns 
gain switched laser pulse. This was a result of plasma rarefaction which lowered both 
plasma temperature and density at late times in the laser—plasma interaction. Plasma 
interferograms taken 15 ns after plasma breakdown showed a low density core ( 

n < 0.1 N. — ) which would have a greatly diminished X-ray emission even with 
continued laser heating and sustained temperature. The low temperature, higher density 
plasma periphery would also be a weak X-ray emitter with the overall consequence of a 


short lived X-ray signal as compared to the laser pulse. 
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Figure 4.4 Electron Temperature for a Gain Switched Laser Irradiated Plasma 
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4.4 Preliminary Observations of Brillouin Scattering 

The first and easiest characterization of the stimulated Brillouin scattering was 
through monitoring the backscattered red-shifted CO, laser light. An 8% sodium chloride 
beamsplitter in the incident laser path collected directly backscattered light (to within =5 
mrad). Earlier laser plasma studies on the polar distribution of CO, laser light scattered 
from this target (Offenberger [53]), determined that during the first 10 nanoseconds 
after gas breakdown, only backscattering occurred. Late time forward refraction of the 
incident light occurred but was unrelated to the Brillouin scattering Besides ray 
re-tracing [54] and a gain maximum [16] to explain the strong directionality of the 
backscattered light, polarization and plasma density gradients could be very important in 
this experiment. Brillouin sidescatter amplification is undiminished for Transverse Electric 
(TE) polarized light and is zero for transverse magnetic (TM) polarization [12]. For this 
reason, the maximum instability gain for an unpolarized gain switched laser beam would 
be in the backwards direction. The TE polarized IML laser beams also exhibited strong 
direct backscatter which suggested that other mechanisms were sufficient to control 
Brillouin sidescattering. These could include the finite focal spot size and the 200 - 300 
micron radial density gradient scale length measured from the interferograms which 
would restrict the plasma length available for sidescatter gain. As mentioned in the theory 
chapter, this is due to the phase mismatch from the electromagnetic wave dispersion in 
the inhomogeneous plasma limiting the spatial amplification. In addition, laser beam 
inhomogeneities could act to enhance the transverse density inhomogeneities and further 
inhibit sidescatter. The cumulative effect of all these features to the laser-plasma 
interaction is that direct backscatter of the Brillouin generated light should dominate. 

Preliminary Brillouin scattered light measurements identified the temporal features 
of the instability. Figure 4.5 shows the incident gain-switched CO, laser, laser light 
transmitted through the plasma and Brillouin backscattered light The relative timing 
between these plots is accurate to better than 1 ns. Axial mode beating in the CO, laser 
output resulted in typically ~ 50% modulation of the beam intensity. Prior to gas 
breakdown, the transmitted light followed the incident light but with breakdown the signal 
fell rapidly to noise level as light absorption and scattering occurred. The Brillouin 


scattered light consistently appeared = 4 ns after breakdown (independent of the incident 
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Figure 4.5 Gain Switched Laser Transmission and SBS 


86 


ay. 


i oe ane it 
of of 
“ ; { 
ws 
a ; 


if 
; { \ \ ' 
: = 1. , 
‘ } } 
' i] f _ 
' ’ j j 
\ < 4 f 
wi 
\ Le 
q 
' 
t 
{ 
: 
i 
‘ 
J \ 
i 
: ~ 
H ~. ts D a i 
F —s 2 mpeg oe ee taney ner Venom errr 
A i Be 2 tea i 
A te 
if th 
é 7 ; 
t= —_" = a 2 ae eae ita Sent etl On «or 


- y 


a 


er" le ae at ee tant Serene 
* 


“Baz 


a 
4 


F 
| 
3) ne 


87 


laser power). This delay seems likely to have resulted from a combination of limited 
ionization rate, laser absorption during ionization and plasma electron density and velocity 
smoothing from hydrodynamic flow. The exact contribution from these was not 
separable from a simple measurement of light transmission. 

The Brillouin signal risetime Ol teans) >> V/o.s ie suggested that 
fast temporal growth to convective saturation did not occur. Time resolved ruby laser 
Thomson scattering measurements of the ion mode confirmed this and will be discussed 
later. The instability growth can be moderated by wave energy loss mechanisms such as 
inverse bremsstrahlung absorption of the pump and backscattered light and ion Landau 
damping or ion trapping acting on the ion mode. However, in view of this long risetime 
and the duration of stimulated Brillouin scattering, a steady state model description is 
likely adequate. Figure 4.5 shows an example of a 4 ns duration Brillouin scattered light 
pulse during which time convective saturation can occur. Over many laser shots, the 
backscattered light varied in pulse width from 2 ns to 8 ns duration with an average of 4 
ns. More quantitative discussion will follow in succeeding sections. One final comment on 
the gain-switched laser is that the scattered light modulation was not completely 
correlated with the pump modulation. Brillouin instability gain is proportional to pump 
intensity, but as seen in Figure 4.5, the scattering terminates before the laser pump is 
removed. This spiky behaviour of SBS has been observed by other researchers ([53], 
[55], [56]) even with temporally smooth laser pulses. No single mechanism has been 
positively identified to cause this abrupt termination of the Brillouin instability. 

Similar studies on the IML laser-—plasma interaction showed notable differences to 
that with the gain switched laser. In the first place, Brillouin backscattering was prompt. In 
Figure 4.6, the incident CO, laser, the Brillouin backscattered light and the ruby Thomson 
scattered signals are illustrated. As seen in this figure, the scattered light turns on 
abruptly within 2 ns of the pump arrival. The ruby Thomson scattered signal provided a 
useful indicator of plasma formation and also confirmed the rapid rise of the Brillouin 
instability. The Ruby Thomson scattered signal had a precursor 750 ps ahead of the 
Brillouin ion mode scattering . This source of this precursor is not known although 
enhanced Rayleigh scattering due to a neutral O resonant state could take place [57]. The 


main point is that the Brillouin instability occurred much sooner after breakdown than for 
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Figure 4.6 SBS and Thomson Scattering From an IML Laser Irradiated Plasma 
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the gain switched CO, laser because of the high initial pump intensities. Again, the Brillouin 
instability decreased abruptly and from Figure 4.6, was clearly not related to lack of 
pump intensity. The instability duration was even less than that normally seen with the gain 
switched laser and may have been related to rapid plasma modification in the strong laser 


pump field. 


4.5 Infrared Light Diagnostics 

In the last section, preliminary remarks were made on the Brillouin instability 
without expanding on its characteristics and how they compare to model predictions. 
Here, infrared spectral measurements and CO, light reflectivity determination will be 
reported and analysed. The complex behaviour of stimulated Brillouin scattering in a real 
experiment is difficult to explain through one self-consistent model. However, a 
composite description can satisfactorily account for the principal observations. The 
results will be interpreted through the ion heating models developed in Chapter 2. All 
together, the observations indicate that the convective saturation of the Brillouin 
instability is principally determined from plasma density gradients and ion acoustic wave 


damping on heated ions. 


4.5.1 Backscattering with a Gain Switched CO, Laser 

The first observations of stimulated Brillouin scattering were made with the gain 
switched CO, laser irradiating the gas target. Red shifted backscattered light was 
focussed onto either a Ge:Au detector or a pyroelectric detector with a suitable number 
of mylar sheet attenuators in front of the detector to ensure linear signal response. 
Reflectivity calibrations were made by inserting a sodium chloride window 
retro-reflector in the CO, laser beam path to provide 8% back-refliected light when the 
laser system was fired into an evacuated target. With average focal spot intensities of 6 
x 1032 W/cm?, power reflectivities of R = 0.20 were recorded and at half this 
intensity R = 0.14 | Yet a further reduction of two in pump intensity to 1.5 x 10” 
W/cm? resulted in a drastic reduction in reflectivity to R = 0.04 . Thus a knee in 
the reflectivity as a function of pump intensity occurred at 3 x 10%? W/cm? above which 


R slowly varied. These results are similar to those reported before on this gas 
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target [24]. The values referred to are peak power reflectivity; as a consequence of the 
short scattering time as compared to the pump duration, the energy reflectivity was 
lower by approximately an order of magnitude. 

Further details on the Brillouin instability, which assist in formulating a model, 
came from time integrated infrared spectra of the backscattered light. For these 
measurements, the backscattered light was directed onto the infrared multichannel 
analyser after passing through a grating spectrograph. Figure 4.7 shows three types of 
spectra which were detected. Figure 4.7(a) is a representatative spectrum of red shifted 
backscattered light generated in stimulated Brillouin scattering. The apparent wavelength 
red shift from 10.59 microns was 55 A with a 20 A FWHM. Before 
deconvolving the instrumental function as discussed in section 3.4.1, the apparent width 
was 30 A , indicating that careful handling of the spectral data was required. 

Several interpretations could possibly account for the spectral width though the 
most probable explanation appears to be ion wave damping. Plasma velocity gradients or 
density gradients could also cause spectral broadening, but are much less likely for our 
conditions. The frequency dispersion due to an approximately 100 micron long Brillouin 
active region, as determined through model reflectivity predictions and direct 
measurement through ruby Thomson scattering , determines the possible spectral spread. 
Doppler shifted ion mode and electromagnetic wave dispersion together alter the ion 
mode frequency by dw % (c. + u)dk + kdu where U_ is the plasma axial 
flow speed . Electromagnetic dispersion enters through wavenumber matching in k 
Without the laser depositing energy in the plasma, it would be simple to estimate the 
flow velocity gradients du/dk from the observed density gradients by using the 
pressure balance equations vn uo/2 = aes . However, the presence 
of both electron temperature and electron density gradients in the plasma compound the 
difficulty of directly applying the pressure balance equation. More rigorous computer 
simulations using CASTOR indicate [51] axial velocities =1 x 10° cm/s which are small 
compared to C. . Also, the measured 600 micron density gradient scale length is 
long compared to the 100 micron Brillouin interaction zone, eliminating it as a likely 
candidate to explain the observed backscattered light spectral widths. Thus, ion acoustic 


wave damping remains as the most likely explanation for the large Brillouin instability 
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frequency bandwidth. 
lon wave damping can broaden the resonant frequency spectrum (even for 
Ak/k << 1 ) and indeed the damping coefficient can be estimated from the 
observed backscattered light spectrum, which we calculate now. Although a significant 
contribution of the backscattered light was due to the rising and falling portions of the 
Brillouin signal, application of steady-state formulae allow spectral width estimates. Pump 


depletion can be ignored for low Brillouin reflectivities whereby 


R= e« exp G where G__ is the convective gain in the unstable region, 

ra 2 2 
G = kK bnw Vo. /8V , vena . The ion wave frequency dependence in 
Gis through the ion wave damping Yo /W. . The noise source € is estimated 
from Equation 2.24 to be CO, laser light reflected from the plasma density gradients. This 
approximation yields € = 107° when n/n. = 0.5 and 


the density gradient scale length is 20 microns. Noise source estimates from the 
reflectivity curves [24] give similar values. Typical reflectivities and pump intensity are, R 


Z 
=0.15 and iV ei) 


~ 


= ] respectively. With an 80 micron interaction 
length and the estimated gain G=5, the inferred ion wave damping from the gain 
expression is Y./W. = 0.6 This large damping factor implies strong ion wave 
dissipation. 

The correct kinetic frequency dependent gain calculation incorporating ion Landau 
damping provides a comparison to the measured spectra shown in Figure 4.7. The gain 


factor G written to include frequency dependence is 


¥./W, 


max (y Jug) 


Equation 2.26 is used to calculate the damping factor, mie as a function of the 
backscattered light wavelength. A curve fitting of the calculated backscatter spectra to 
the experimental spectra in Figure 4.7 using the proper susceptibilities for variable 

ie ; T, gives the plasma electron temperature and electron—-to-ion 
temperature ratio. The calculated spectral peaks were sensitive to the electron 
temperature while the widths were sensitive to the temperature ratio. However, the 


spectral shape was not strongly affected by the choice of Ge . This 


ad? iio? ugteciiiae st ries jvigtncdideiy peieisides ent | 
Pantie. s Hevorttis wight aamiahga. by. Tai: On 
ai 6. enicinicG orillet be goign ter 6 At avy new 


anu" astamias tbiw ihvioage wgiy eens + tye rh yoi ro rosso arg 
yaa erin aciviesRite rises}, . yaal 1s Saupe #0 ma : 


» 


fieae eideient Sf 47 SQ AVIA VUrIGoO art Ee oD Basi 8 sal "1 ca 


phaach 52 2 BON amele xt i 4 A re heanaaeeeat ces 
HHIT.2INsibaD yHenst eee early tI 77 bsiasitiah a nat (2-01 oF 85 

brig | a.0 = Who Tay "54 at. oe abbey ¢ e are w) 

wit (mor. 2etanii2s =i wee Asi 2h6 rien OS 2 ethic aison rowing veo 

| ghtva brs esisivirselton isticyt apie te en 

rneietel domain’ Os nd ani evinssdast® oh ‘: They ty us > - 
ep. art: mont. onigrisb. avew ete a ane sish a0 ng ane pe ag : a 

avew lol Baad esi ic ieoh at vests sie hod ol ihe air 


keg. Be i ra —s 


oo ‘st ian esa niamto wtf a 


5. aqugi MY quota onsege oe neaaqmen « 2ebiverg | | 
| is “ponoupent seat orca Oe 


pideiisy 70% coiliinannaue- saat seo Geta oa unt ri “evaega mene 
Noinet=navisela ins sure pageret nigrpels arratig oett sovig 
navosls Brt ot SBvitians2 sev, eg, shied Cera OID aft piger 
ant evawor oiler wren fe ote: vidas sew aritinw ade : 
dal xen? a: a | | 


i” 7 it m =) a 
om. ) — 


93 


weak Gay dependence increases the confidence in the curve fitting estimates 
of ss and ZT /Ti from the backscatter spectra. The dashed line in Figure 
4.7\a) is a curve fitting with 3 =130 eV and ZT /T; = 1.0. This value for 

T is in good agreement with the value of 120 eV determined from X-ray 


measurements. The large ion temperatures 0(ZT,) , however, are not expected from 
electron-ion collisional equilibration alone, further justifying an ion heating model. 

We now address the question of ion heating. At low Brillouin reflectivities, heating 
estimates from the energy deposited to the ions in the Brillouin interaction show that 
there may not be sufficient time to heat all the ions. The fastest ion heating rate for a 


stationary plasma would be 


dT, /dt = (w./w,)(RI/n.L) 


where w/w, is the Manley Rowe energy fraction transferred to the ions 
in the Brillouin instability, R is the power reflectivity, I is the pump intensity, 
ns is the ion density and L is the interaction length. For 20% Brillouin 


reflectivities and 1 x 10!5 W/cm? pump intensities, the =600 eV/ns ion heating rate 
calculated from the above can account for the high ion temperatures. However, for 5% 
reflectivities and 3 x 10'? W/cm’, the ion heating rate is too slow. Under these low 
reflectivity conditions, it seems more probable that some ion fraction is trapped and 
heated to eee O(ZT.) . Indeed, simulations would suggest that this 
occurs [11]. The broadened gain bandwidth due to ion Landau damping on the fraction of 
heated ions n,/n , would be the same as if all the ions were heated. The 
damping factor, however, is diminished by n,/n and the curve fitted ion 
temperature estimate is that of the trapped and heated ions. Any low temperature ions in 
the plasma would not broaden the Brillouin backscatter resonance. 

The spectral data are not always so easily interpreted as seen in Figure 4.7(b) and 
47(c). Both of these spectra show significant specular reflection as well as the clearly 
identified Brillouin backscattered component. In Figure 4.7(b), the narrow ion wave 
feature is curve fitted as above only now the electron temperature and electron-to-ion 


temperature ratio are sis = 90 eV, ZT /T; = 2.5 . The smaller ion Landau 


sxupit nh ani! penta ant gitaaae » ogame sine mW GS | 
Ot Sev aT De. yt "hee vin has Br a 
yer“R mort wenimaiet Ve ORT 1o suliey atit 4 ui 
seit betreane for ana revewwdel nel, T7009 i | 
lehoen gadtken cal me BAN ee eT malig | 

er asiviteet ie} Telos: & dae th sgniiaer eh ie neta casi at 
igi Worle aoitcsreie migah) ss sre Gt 4 at ch agiongen yorene. eet mart 
; Wa? sienigmisgor! Adi f2efagt anT (eng ott " seer of ee yaigid tue vl 
| ed Sivow a 


i ; ‘ Pe. 


: 1p ONEADK Ae oe 


| ahi a oo a , 
wt oF te retenad norset ve ‘ane owe tonne 4 a 
rianated oruQcadr et: Doty soeites Tenia ‘sit ai a etifielecen 
| " a. 
Wold 20S oF “dtpasl ‘nomoe et “ent ae _ ei lng eat A 


~ 4 
oo a sen 
sigy onitee dt sriive-GOgs ad aaeggi am yon re a bre 2 ad 
9G not leVeWwOoH sate ied weirs | 
WO|l Beart wool Wwole’ oor si 
so) atl +y a" t *) 4 why § 44 = 


ir? jeri teenpwWe oMOW 


Ouse Sitio See Wet bog ok he 3 


i” 
2 
co 
a 
‘G 
2 
mM 
=o 
oe 
@ 
st - 
—_- 
~ a 
cm 
~ 
a 
is 


noi beth? aveia ‘srt ola 
mi Arvo suite 149 eriet wo vA eat ) 
eanstigaat 


ors lit Ss a tEQIATH saan 26 eet 


yield sot .ae lisw as noiragite9 


a mews os <n omg a ae 


avaw AiG) wor var, art Caatly. 


vaBnegd 101 rallatria or: 


a | 


94 


damping is attributed to decreased Brillouin reflectivity resulting from a lowered pump 
intensity. In some laser shots, the time integrated backscattered light had a large specular 
reflection coefficient. Figure 4.7(b) and 4.7(c) are two examples of this. Overdense 
plasma could form in the supersonic gas target during the laser—plasma interaction and 
reflect incoming CO, laser light. This specular scattering from overdense plasma and the 
Brillouin scattering do not necessarily occur simultaneously or at the same power levels. 
Unfortunately, the temporal features of the specular light component are not known 
from these time integrated spectral measurements. Figure 4.7(c) also has a small amount 
of 20 A blue shifted light due to Doppler shifted reflection from a moving 
plasma boundary. Overall, the backscatter spectra in Figures 4.7(b) and 4.7(c) occurred 
infrequently and the majority of the backscatter spectra were well behaved Brillouin 
spectra as illustrated in Figure 4.7/(a). 

Although the reflectivity and spectra of the Brillouin backscattered light could be 
accounted for through ion trapping, heating and Landau damping, the effective interaction 
length, remains to be explained. For a Gaussian laser beam, a diffraction limited focus has 
a half intensity waist A/t x (focal length/beam radius)? = 240 microns. This 
calculated beam waist dimension is within a factor of 3 of L calculated from the Brillouin 
gain and that determined from ruby Thomson scattering off the ion mode, but it is 
misleading. As explained in Chapter 3, the focussing parabola’s aberrations extended the 
effective beam waist to 1 mm so that the interaction region was not determined soley 
from the pump intensity. 

Density gradients in the plasma can potentially limit the resonance zone for the 
instability. Inspection of Equation 2.18 shows that when H(x) »v 0(1) , sufficient 
dephasing of the coupled waves reduces the instability gain to zero. For an assumed 
linearly varying wavenumber mismatch resulting from electromagnetic wave dispersion, 


the estimated unstable region is then: 


eer nil a 
n C 
ec ee 
ky n/n. 
where LY, is the density gradient scale length. This limit on the unstable region 


occurs because the ion acoustic waves can propagate into plasma where they no longer 
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satisfy the local 2k, matching with the electromagnetic wave pump. From the 
interferometric measurements showing LY =600 microns and n/n. =0.5, we 
estimate L =50 microns. This interaction length is less than the 80 microns 


determined from reflectivity and ruby Thomson scattering measurements. It is likely that 
during the interaction, a number of ion modes are excited whose resonant points are 
axially separated. Once the plasma density gradients fix L, the heavily damped Brillouin 
instability convectively grows as described by the ion heating or ion trapping and heating 
models. In summary, ion heating and density inhomogeneity appear adequate to account 


for the spectral and gain length characteristics observed experimentally. 


4.5.2 IML Laser Backscatter Spectra 

Backscattered light observations were made with the IML CO, laser as a pump 
source. The large pump intensity range was sufficient to give a wide range of Brillouin 
backscatter reflectivities. Figure 4.8 shows the reflectivity as a function of laser intensity 
over the intensity range 1 x 10% W/em? < | < _— 1.7 x 1033 W/cem?. The nonlinear 
reflectivity trend can nominally be explained through ion trapping and heating. Application 
of the steady state model of section 2.3.2 qualitatively predicts a reflectivity minimum 
because ion Landau damping does not increase significantly once all the ions are trapped. 
The steady state ion heating model also exhibited a similar behaviour but because of the 
long time required to heat all the ions ( Ge. ZKT Nn. Lu /RIw, ), it was not 
likely to occur, at least for short interaction times and low reflectivities. The dashed line 
in Figure 4.8 was computed from the ion trapping and heating Equation 2.32. The best 
fitting plasma and interaction parameters for these results are, € =10%, 

7 le = 15, n/n. =0.3 ¥o/M. =0.3. The large electron-to-ion 
temperature ratio is that of the plasma before ion heating from the Brillouin instability. 
The fitted curve shows higher than measured reflection. This could be in part due to the 
assumption that only trapped and heated ions damp the ion wave. The ion energy 
equipartition time is 400 ps and so some heating of all the ions can take place during the 
laser—plasma interaction. Increased ion Landau damping of the Brillouin ion mode would 
result and the measured SBS reflectivity would be lower than the theoretical estimate. 


Once the reflectivity begins to increase again at I= 3 x 10’? W/cm’, almost all of the ions 
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Figure 4.8 SBS Reflectivity from an IML Laser Irradiated Plasma 
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are trapped and ion thermal equilibration is not important. Above this pump intensity, the 
Brillouin ion mode density fluctuations approach unity and nonlinear effects such as wave 
breaking (Forslund et al. [11]) can take place. Nevertheless, reasonable agreement is 
shown between experiment and model calculations based on ion heating. to further limit 
the Brillouin instability. 

The time integrated spectra of the Brillouin backscattered light was still broad as 
expected for ion Landau damping on the high temperature ion fraction. Figure 4.9 shows 
a series of infrared spectra recorded for conditions OrierReenO 3-539 _ At 
these intensities, the pump level determines the instability strength and very large 
reflectivities were observed. Specularly reflected infrared light was not detected 
because the short laser pulse did not have time to ionize the gas target to critical 
densities. 

Theoretical spectral curves were generated by the same method as those for the 
gain switched laser backscatter spectra. In all three theoretical spectra plotted in Figure 
4.9, the electron-to-ion temperature, ZT O/T; , was selected to be 1. The ion 
temperatures corresponding to these theoretical curves are 120 eV, 240 eV and 300 eV 
in Figures 4.9(a), 4.9(b) and 4.9(c) respectively. Together with the measured 150 + 20 
eV electron temperature, these ion temperatures indicate a low average ionization. If 

T. . 0(ZT,) from ion trapping and heating, the average ionization, Z, is less than 
2. However, interferometric measurements of the electron density in the IML laser 
irradiated plasma indicated Z > 3 . This contradiction in Z suggests that the ion 
temperature is less than ZT, , although a decreased T, would not account for 
the backscatter spectral width. Other laser-plasma effects such as Doppler shifting from 
plasma motion and ion mode frequency pulling due to the intense laser pump (see section 
2.2) could have contributed to altering the time integrated spectra from theoretical 
estimates and making their exact interpretation difficult. Also, unlike the long pulse 
spectra, the short pulse spectra exhibit tails to the longer wavelengths, suggesting that 
the differences may be related to time dependent conditions for the short pulse case. 

In summary, these backscatter spectra demonstrate again the elevated ion 
temperatures and strong ion wave damping in the laser-plasma_ interaction. The 


quantitative agreement with our steady-state model is not as complete as time dependent 


ary iienpin que aid b avods fei de 


vsw ae daue ayaetle eae gg we os idinendea nib bangs. 
fim) Wari wT Ot aniread ri9 46 oe snpitatiaisa at eet bi 5 InomiwoRe / cn to 8 " aN ned rwor 
oo - hs ai ne : it : i’ i 

, ZW trp erectenee pigeitine) eet 1 siroege hinvigstl wr ot 

ote e+ Siu A sex AO} seiner? dnt seth arty rao. griginiks vabrad nol iat B 

= Ae 1g. | _anigitiotws 10? deb eoRT wapege tr tod 

‘Ne we ith orl eéhicinerey aval tyrud ate 

caw | Siok pereyin Setoertan yiralIege baviwedo oyew- 


Qual 28g) a , ak ot nt ati eyed ton bib ealug cal note ag 


te 


‘A1armn. 3 mee ort re bstprerag sow seirunieiabel ot 
a et sq (sal ie toa ent te, fit site vategae ined Vent bn 

parcels eth ae ee soroi-nowoa 
O+S ve OST me ea\qus) lasircibit: saert'ct LIES Bel 
arn att ‘top s yomnee ose 
be msec cee wear swigreginat rh’ 
srit geal ai Syrtoiant preg Bak gniggey envio >, sisi * ce 

gost JIM get om wiieees ne ocd a oy “ ’ cain nl 


Notaesinor. epereve weal se _stealbn 


rae) t 


1977 OA ittigas sigqoG 2 oe 2 {v1 
1olte2 382) ermwa, ieash eanatan ant or Eateasien Sbom Mot bre roitom eral, 
jacitergsnt, Mott sic oe borergginl smi ei paar at vetudinings on h ss ae 
a2lg gach Bn iin oe Borin ‘testi: rea saath soaKe inst am ons aaitaedlical iu 
ise? unitéeppue aegrisigvaw rego! edt < \ flat iin, ered ‘sehug ode eit aitoege rs 


Ts) 


S252 sziug morte art Bro sndigibn@S ropragin ani baraler 34  anlasacatiints At ine 
_20rt yneaneme a ee aan 


noi Pejayaie sri nisgs sie vencmab, sutosge, 18 
‘oe 


“a v0, Bend hia mnt 


aif “noithansini siwailg- aaa ait. Ai _ phic 


nebneqah anit 2B etalon: 26 tan gue eit 2. tC ae ‘i 
J i" ‘ 74 oe * 


j mi. Fy y i Se 


98 


l= 1 x10'3 wyem 2 
R=0-53 


Jot ties ly bei ~~ 
do ho + 50A 


[= 5x10l2 w/em2 
R=0:23 


(C) l= 5x 10!2w/em2 
f R=0-10 


wi hot SOA 


Figure 4.9 SBS Spectra for an IML Laser 
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effects are likely of some importance for short pulse irradiation. 


4.6 Stimulated Brillouin Scattering with a Multiline CO, Laser 

Stimulated Brillouin backscattered light levels are a direct measure of the 
instability’s strength. In view of its practical importance in laser fusion, a series of 
experiments were performed in an attempt to decrease Brillouin gain and consequently 
reduce the amplitude of the backscattered light. As described in Chapter 3, both the gain 
switched CO, laser and the injection mode-locked laser could be induced to oscillate on 
two or more laser lines simultaneously. When the laser frequency separation Aw was 
WZ 


Vos py! 2e : 


our experimental results verified the theoretical predictions of reduced SBS with 


much greater than the homogeneous growth rate y, = (w/w) 


increased number of lines as expected. More complex behaviour was observed when 
Aw Xe for which backscattering from the various laser pump lines were 
not decoupled. 

Brillouin backscatter reflectivities with the gain switched CO, laser operated both 
single wavelength and multiline are plotted in Figure 4.10. Average single wavelength 
(10.59 microns) power reflectivities were seven to ten times higher than that observed 
for the same total pump intensity equally distributed in 10.59 micron and 10.55 micron 
laser lines. The residual 2% reflectivity for more than one line did not change with three 
or four equal intensity lines and, moreover, showed longer temporal features than the 
Brillouin backscatter light. These observations indicated that this remaining light was 
specularly reflected from a critical density layer or steep density gradient somewhere in 
the plasma. This behaviour is consistent with the spectral measurements mentioned earlier 
as being unrelated to the Brillouin instability. The Brillouin gain reduction was greater than 
the algebraic reduction expected from the ion heating model. In these results the 
reflectivity reduction is closer to an exponential reduction R ~ e exp(G/N) 
when modifications due to ion heating are not incorporated. For G =5, two lines 
result in a 12 times reduction of the reflected light. From this, it is very likely that the ion 
heating model is no longer valid at these reduced reflectivities since the ion heating rate 


is correspondingly low. 
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Figure 4.10 SBS Refiectivity for a Multiline Laser 
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In yet other experiments with a multiline pump laser satisfying Aw % a 
the SBS reduction was considerably less than above. For roughly equal line intensities at 
10.289 microns and 10.274 microns for a total irradiance of 5 x 10}? W/cm?, the SBS 
reflectivity, shared equally between the two lines, was reduced only 30% from single line 
lasing at 10.289 micron. The small effect of the two line irradiation on SBS in this case is 
attributed to the lack of wave decoupling when sy Yo 

Further evidence of anomalous behaviour when Aw es was seen when 
one of the laser lines was much weaker than the other. Figure 4.11(a) shows the time 
integrated spectral content obtained with the infrared multichannel analyser. This plot 
shows that the SBS reflectivity of the weaker line approached unity whereas the 
dominant line maintained =9% reflectivity. For this particular shot, the line intensity ratio 
was 1:11 and the SBS ratio was 1.2:1, indicating that, within experimental error, the weak 
line was entirely Brillouin reflected. This suggested that the strong pump drove ion 
fluctuations from which the weak pump could scatter. When a third pump line at 10.26 
micron was introduced, the presence of the very weak 10.274 micron line led to only 
moderate reduction in SBS reflectivity. In Figure 4.11(b), the laser intensity was nearly 
equally divided into two lines where Aw was large enough to expect significant 
reduction in reflectivity. However, a weak 10.274 micron pump was strongly reflected as 
it apparently coupled to the ion modes driven by the 10.289 micron pump, consequently 
there was no reduction in overall reflectivity. In this example, the 10.26 micron line did 
not produce Brillouin scattering in spite of it being much more intense than the 10.274 
micron line. This clearly shows the complex behaviour of SBS which can result with 
multiline irradiance when the line separation Aw does not significantly. exceed the 
homogeneous growth rate. 

In order to determine whether any significant temporal differences in SBS existed 
for multiline irradiation, the gain-switched CO, laser was operated at 10.28 microns and 
8.55 microns. For these widely separated frequencies, narrowband filters and two Ge:Au 
detectors were used to monitor the individual laser lines or their associated Brillouin 
backscatter. As expected for these two widely separated lines ( 


Aw = 1.4 x 10! 3 aa gah ), the reflectivity was greatly reduced. Within 


the 2 ns instrumental resolution, both lines exhibited identical characteristics in input and 


FOF 


i's aoirenerri ori lines ghiguen 197 evoda rari axel: ie abianes save nowsuber® 
See ort “meivy Ol x & te ganateets lato? @ wat aroha BSL brs enorain & 
ori gignie most 0L ying besubiey eew..corit ows set reainied viiaups barerte 
+i gee avit oi 2B2 me nodsibert e4-oud ait te ‘selte line sa? novo 88.0% te 


meg ¢ 


org fuani mi 2diteineianiss < eseitng sont A mong aree seit al 8 


ae nn ee ee euotsenane' Fo sonsbive rer oa 


Ay ocmeviens fennatiodigm GCaxaiirnh oe ar iw bernie InistrioS iesoage t 
saatw ying pertonc ge srt. w sav oh % vtvitaathan ea ort wre 
viene) ard) eft Joe. sahioined-sint 164 iivitgahien OF banigsri . 


9 282m egonare thd ie roaner trsci Piewgla’ ra wemeriw eniecienels oF vabioinl 


) owt bos 216i? bredivo rin 2erorbene " oapemgee yiebhn avert 74 anaroin aae ul 
HiNG. bstsiae ceeal Ne 20" sani Bes): jeubiwibai get NONIoN of beau anew ar0t00teb 
earl 087618 mae viebiw owt cert ) -” poaess 2h retteoetoad 


¢ i na on eR wrnig arity e Athi: . 


eo 


{5 


_ 2 vie 


sr * os penleuaady paw to oat salt oP b 


pore jell {[ & ewert cere - eee ont story eB ‘29nil neta! on to 


[he a 


smiegkes rirtiue Jeet QritRotont , 1 aaa widen 282 act 

“Tt: OMorre art? gaat be aeague awit betwen muta ol 

Mud Cael 6 rp yy wirese tiles ara Aaa “eo nich mort oie 

TS.0! dgew wey ent Fo vOngee ty ovis _beobo oni 

Henoin ageel ari. Sal al yt vitagitay ‘eae A oitoubes one 

JeQXS OF. Modorer opie RB ache ‘mage ean awd ‘etn nedivit 
site. vine 2 eeu (amtue:n ity + ETS aera | siag.2 “eweweH Pivitgarter (i ad 

sénea \qmuq colsinn BBS Ol) sree svn seer Ae ext et belquos (unaregge' 

orovh SS.0t’ HF lelumexe iat | nfs _omigattin Neraxe Hi neitoybe” on enw | 

y ned? eanerAl atom igo rn 1 edie ti: to ory 1) jormretiaca rigor eoubexg 0 

Hon S82 to ruojvedsg eaierned att @worta yrealo aitT cenit novoin 

; 2806! Agia rents "Bae: anil emt pervs aoreing yy enilisiarn = 
i ate cnwory suosn@Romod . 

sm BSOt ts balsas aay rou tA fray baripsiwe-nisg ert Noms Tt enilithur ot / 


= 


Ret = Suber (Bea saw vlivooaliss art if a ei ba Eloy xd, | #) tad 


(a) 


Incident 


SBS 


Incident 


SBS 


102 


10.260 um 10.274 um 10.289 um 

cs) 
NN 

\ 

1 & 

6 1 g 

i. 

Pe Noise Limited ‘age 

\ \ 


10.260 «wm 10.274 um 10.289 um 


Figure 4.11 Multiline CO, Laser Backscatter Spectra 
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SBS reflection characteristics, thus precluding scattering from the different lines at 
different times during the laser—plasma interaction. 

Finally, SBS reflectivity measurements were made with the injection mode-locked 
CO, laser to verify similar SBS reduction for short laser pulses. Unfortunately, additional 
laser cavity losses from the intracavity SF, absorption cell lowered the available focussed 
intensity to a region where the SBS reflectivity was already low. A maximum intensity of 
2.8 x 10! W/cm? was achieved from which a single line (10.289 micron) reflectivity of 
20 + .3% was found. With the same total intensity in two lines (10.289 micron and 
10.274 micron), the reflectivity was reduced to 1.2 + .6% and the three line (10.289, 
10.274 and 10.260 micron) reflectivity was decreased to0.5 + .1%. As with the gain 
switched laser results, it was likely that the reduction of SBS reflectivity was somewhat 
moderated because of the small frequency separation ( 


Mieotex1 0 sn! y =10"¢ elie y 


consistent with the expected reduction of Brillouin backscatter reflectivity for a multiline 


Nevertheless, these observations were again 


laser pump. 

In summary, a dramatic reduction in Brillouin reflectivity has been observed for 
long pulse multiline laser irradiation of a plasma. This is the first experimental 
demonstration of this very effective method for virtually eliminating the Brillouin process 


in a laser—plasma interaction. 


4.7 Ruby Thomson Scattering Measurements of SBS lon Modes 

The Brillouin instability excites both electromagnetic and electrostatic waves in a 
plasma. The backscattered light has been analysed showing the essential-details of the 
instability, that is, ion trapping and heating along with convective growth and plasma 
inhomogeneity moderating the growth rates and saturation levels. Here, the ion acoustic 
wave is examined directly to further clarify the behaviour of the instability. The diagnostic 
method employed was ruby laser Thomson scattering, which has been described in 
Chapter 3. Scattered ruby laser light was detected with either an Optical Multichannel 
Analyser for time integrated spectral measurements, a photomultiplier for average 
scattered light measurements and absolute fluctuation level determination, or a high 


speed streak camera for temporal details of the ion wave k spectrum and its spatial. 
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distribution in the plasma. 


4.7.1 Determination of Absolute Level of lon Fluctuations 

Probably the most convincing evidence for the existence of a heavily damped 
convective instability was obtained from the photomultiplier scattered light 
measurements. By determining the ion mode fluctuation level for known pump intensity 
and backscatter reflectivity, an effective damping constant was found which directly 
demonstrated strong damping. The observed simple functional dependence of the ion 
fluctuation level on reflectivity, én/n « VR , Was the signature of heavy ion 
wave damping. The ion fluctuation level ¢n/n was Calculated from the Thomson 
scattered light power. 

Calibration of the detector response was made with a tungsten filament 
blackbody source positioned at the CO, laser focus. Figure 4.12 shows the two lens 
scattered light collection used in the photomultiplier detection system. The tungsten 
ribbon filament was placed at the parabola’s focus and oriented normal to the collection 
axis. The blackbody emitting area was determined with an iris located between the two 
lenses. Once the filament temperature T was measured with an optical pyrometer, the 


blackbody source light reaching the detector could be calculated from 


P = I(A,T)A da da 


where |, the blackbody emissivity is 


1 Od 9 ac (2he“/a”) (2) Lexp(co/ aT - la | erg/A S cm? sy 
Cy = 1.438 °K cm 
e = 0.43 for Tungsten at 6943 A 


and A is the effective emitting area, dQ the collection solid angle and di was 


the detection bandwidth. Here dA = 4A was determined by a narrow 
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Figure 4.12 Blackbody Calibration of Thomson Scattering Photomultiplier Detector 
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bandpass filter located in front of the photomultiplier. This calibration together with the 
known attenuation introduced to obtain linear response to scattered light signals allowed 
absolute Thomson scattered power levels to be calculated. 


The scattering form factor S(k) which is related to the ion fluctuation level 


2 
| sn| could be deduced from the Thomson scattered power 
p In Vo. da 
S(k) = P scattered! Pee 
where P is the scattered light power, No the plasma electron 


scattered 
density, V the scattering volume estimated from the intersection of the ruby probe 


beam and the Brillouin active region, I the probe beam intensity, De the Thomson 
scattering cross section (8.8 x 10-%* cm? ). The scattering solid angle d2* was not 
necessarily the same as dQ, the collection optics solid angle. As will be discussed 
later, the k spectral width of ion fluctuations was measured to be <Ak/k = 0.20 
which translated into a small angular spread of =0.6 degrees in the scattered light. 
Consequently, the effective Thomson scattering solid angle at the first collection lens 
was less than the available collection solid angle. 

Once S(k) was calculated, the fluctuation level, én/n could be 


determined from the Thomson scattering theory result [58]: 


sn (1 32 i 
—~ = {—| |—| {S(k)k™ dk 
Me OT en 
e 
if 3S(k) is assumed to vary slowly over the wavenumbers of interest, then defining 
(25 A)k, a Keene d)k, where A is a measure of the k spectral 
width, we have: 
én2—s «164 
— = —— S(k) 
2 3 
n x A 
where i is the wavelength of the CO, laser light inducing the collective ion 


fluctuation. 
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The experimental results for an average gain switched CO, laser intensity of 3 x 
10}? W/cm? are plotted in Figure 4.13 as (én/n)2 versus backscattered light 
reflectivity, R. The ion fluctuation level is seen to be linearly related to the backscatter 
reflectivity as expected for a strongly damped convective instability. In steady state, 
Drake et al [9] have shown have shown that 6én/n is related to the 


ponderomotive force by 


én = : 
pera Ree tata pe gene 
2 =O -0- 
+ x. 
n 8M (I ie ie x2) 


which reduces ( as shown for the coupled mode equation calculations of Chapter 2) to 


ee: Yos!2 &s | 


n 2 ey. 


From these equations and a curve fit to the results in Figure 4.13, the damping 
coefficient, A , was estimated to be 0.3. The solid line in Figure 4.13 is the fitted 
curve. This large damping coefficient is still less than Y,/W, ~~ 0.6 which was 
estimated from the Brillouin backscatter spectra. It is difficult to isolate the reasons for 
the different damping coefficients although it must be realized that all the models used to 
obtain these figures are approximate. Possibly the lower reflectivity ( R ~ 10%) in 
these results account for less heating and a lower Ye . In any case, both these 
damping coefficients are the result of Brillouin induced ion heating and ion Landau 


damping of the ion acoustic waves. 


4.7.2 Spatially Resolved Thomson Scattering 

The remaining Thomson scattering experiments used the high speed streak 
camera for light detection. The streak camera not only had the advantage of providing up 
to 10 ps temporal! resolution, but the multichannel detection capability allowed streaking 
of either spatially imaged, wavenumber resolved or wavelength resolved Thomson 
scattered light. Each of these measurements will be discussed here. In the first two 


arrangements, a two lens telescope was used to relay Thomson scattered light to the 
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Figure 4.13 SBS lon Fluctuation Levels Determined from Thomson Scattering 


108 


SCHAHA 


4 


i 


i 


cpo 


: A i A ; 


109 


streak camera slit. For spectral measurements, the light was dispersed through an E1 
monochromator located at the focus of the second lens and its output was imaged onto 
the streak camera slit with a third lens. 

The length of the stimulated Brillouin scattering zone is an important parameter in 
describing the instability. The infrared backscattered light levels have been analysed from 
which an interaction length is inferred. A rough estimate of L from the Bragg reflectivity 


relation is 


(L/a)® = R/L(n/2)°(n/n,)*(8n/n)*] 


Here the ion fluctuation 6én/n is not self consistently calculated from the coupled 
mode equations, but rather is determined experimentally. With the gain switched CO, 
laser parameters, this equation implied 1 esa for (sn/n)? = 0.02 
which is much smaller than both the CO, laser focal depth ( 500 micron) and the 
density gradient scale length ( Ly = 600 micron). 

Ruby Thomson scattering from the Brillouin driven ion modes allowed direct 
measurement of the stimulated Brillouin interaction length as well as the fluctuation 
distribution. The f/10 focussing optics’ 6 micron resolution limit did not interfere 
significantly with the spatial imaging. For both the gain switched CO, laser and the 
injection mode locked laser, the average interaction region FWHM was measured to be 
80 + 20 microns although a range from 50 microns to 160 microns was seen for 
different laser shots. Figure 4.14 illustrates a streak camera plot (relative scattering levels 
only), of SBS ion fluctuations in a gain switched CO, laser—plasma interaction. Thus the 


-1/2 


directly measured value is greater than either the inhornmogeneous L v k' or the 
Bragg L_ but is less than the optical depth. This gain switched laser result was 
particularly interesting in that it had two intensity peaks 1 ns after the most intense phase 
of the interaction. This is supportive of density gradients or at least axial density 
nonuniformity determining the instability zone. The spatial density fluctuation distribution 


at the time of most intense scattering is qualitatively similar to the damping calculation of 


section 2.3.2. It must be remembered that the wavenumber dephasing from a density 
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Figure 4.14 Time Evolution of the Spatial Distribution of SBS lon Fluctuations (Gain 


switched Laser) 
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gradient alters the coupling coefficients in the coupled mode equations. Two examples 
of the ion fluctuation spatial distribution for an IML laser pump are illustrated in Figure 
4.15. The two plots show the observed variation in the interaction length. The most 
significant difference between these results and those obtained for the gain switched 
laser is the shorter time duration of the ion fluctuations. 

The 80 micron average spatial dimension of the scattering region corresponds to 

c= 5 in the steady state, WKB spatial solution of Equation 2.20. 
For the laser and plasma conditions, the corresponding normalized damping, 8 is =10 
from which one infers Held = ] . This anomalously large damping figure is 
probably due to inexact determination of the laser pump intensity in the interaction zone. 
In any case, the large ion wave damping inferred confirms the convective gain model 
concluded from SBS reflectivity and backscattered light spectra as well as time 


integrated Thomson scattering results. 


4.7.3 SBS lon Acoustic Fluctuation Wavenumber 

The ion fluctuation calculation from the photomultiplier results required an 
estimate of the wavenumber distribution. This figure was obtained from the streak 
camera measurements of the wavenumber resolved Thomson scattered light. The streak 
Camera entrance was moved away from the second lens focus of the collection 
telescope used in the spatial imaging experiment. As discussed in Chapter 3, the light 
detected across the streak camera slit was proportional to the Thomson scattering angle. 
= 2k 


From Thomson scattering theory, k sino/2 or for the small 


k 


ion ruby 


scattering angle © probing Brillouin fluctuations, _k . Hence the ion 


ssa 0 
ion ruby 
wavenumber on was immediately calculated from the scattering angle measurement. A 


typical wavenumber spectrum is plotted in Figure 4.16. The wavenumber width, 


Ak/k. 4, = 0.15 in this case is slightly lower than the average, 


Ak/K. oy, 


calculated for k corresponding to the vacuum CO, laser wavelength because of the 


0.20 . The scattering angle was reduced from the 7.5 degrees 


uN 


plasma refractive index, i.e. 


k= 205" - n/n. /¢ 
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Figure 4.15 Time Evolution of the Spatial Distribution of SBS lon Fluctuations (IML Laser) — 
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Figure 4.16 Wavenumber Spectrum of SBS lon Fluctuations Determined from Thomson 
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A more careful calibration of the absolute scattering angle would evidently allow the 
plasma density in the Brillouin instability region to be calculated. 
A significant observation also made was that no subharmonics ( Ky ) or 
harmonics ( 4k, ) were detected with the pump levels used in this experiment, 
I< 6x 10? W/cm? even though saturated Brillouin backscatter levels of 10 to 
20% occurred. Either of these harmonic fluctuations, if excited by the Brillouin decay 
products could act to saturate the instability [59] [60]. This absence of harmonic ion 
waves contrasts sharply with experiments of Walsh and Baldis [61] where 4k, was 
seen even for 6 to 10% saturated backscatter reflectivity with a short pulse CO, laser 
pump. Both lower pump intensities and higher ion wave damping can explain the 
difference in results, though the stronger damping is a more likely cause. This result again 
confirms the heavy ion wave damping of the present experiments. The absence of 
subharmonic spectra rules out secondary ion wave decay as a saturation mechanism in 
this case as well . This too was consistent with the large ion wave damping since decay 


requires én/n > Y./¥. 


4.7.4 Temporally Resolved SBS Ion Fluctuations 

The 10 ps temporal resolution of the Hammamatsu streak camera allowed 
accurate determination of the Brillouin instability growth as well as decay along with the 
observation of high frequency fine structure, if present. Figure 4.17 shows typical streak 
camera results of the spatially imaged Thomson scattered light with a gain switched CO, 
laser pump. These figures illustrate the monotonic growth of the ion acoustic fluctuation 
level which over many shots showed growth periods of 150 to 700 ps with an average 
of 500 ps. The instability decayed in similar periods. The absence of exponential growth 
in time or rapid modulation in the Thomson scattered light was clear evidence for strong 
ion wave damping appearing well ahead of the peak ion fluctuation level. The scattering 
risetime was less with the IML CO, laser pump (average risetime =400 ps) and again did 
not show rapid modulation. This temporal behaviour was confirmed over a detection 
dynamic range of better than 10,000, obtained by using different neutral density filters 


to attenuate the Thomson scattered light and by changing the camera streak speed. 
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Figure 4.17 SBS lon Fluctuation Growth (Gain Switched Laser) 
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The homogeneous growth rate Y gives a risetime of 1 ps, or for heavy 


O 
damping, simple convective gain calculations give a risetime estimate of GL/c =12 
ps which are far too rapid for the observed signal growths. The more detailed 
convective growth rate calculations from the coupled wave equations in section 2.3.3 
also do not seem to apply. to the above observations except that the lack of any 
oscillation in the ion fluctuation level is qualitatively consistent with the strongly damped 
wave Calculation. 

On the other hand, ion trapping and heating is a wave energy sink that can account 
for the slow rise to saturation. As the ion wave amplitude grows, energy is expended in 
trapping and heating ions. From balancing the fraction of the incident pump energy 
imparted to the ion acoustic wave and thereby to the trapped ion energy, a growth rate 
can be calculated. In this model, the instability convectively saturates to levels determined 
by ion Landau damping and/or plasma inhomogeneity. Now, Kartunnen and Heikkinen [62] 
previously demonstrated, by solving the coupled mode equations including wave damping 
through trapping, that the saturated instability levels were unaffected by trapping, even 
though the growth rates were . Their formulation did not include other wave damping or 
account for wavenumber mismatches in inhomogeneous plasma. Conversely, Rosenbluth 
et al [7] who solved the equations including both plasma inhomogeneity and wave 
damping but ignored ion trapping, did not find a reduced growth rate to saturation of the 
instability. Furthermore both of these models did not allow for increased ion Landau 
damping on the trapped and heated ions. 

Such a model is proposed here. A steady state solution to this problem was 
derived earlier in Chapter 2. A rough estimate of the Brillouin instability temporal growth 
was made from the energy balance which ignored ion convection and trapping dynamics. 
Three equations were used to obtain an expression relating the trapping velocity width 


and the Brillouin backscatter reflectivity. These were 
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all of which have been introduced earlier. Eliminating ¢ and 6n/n_ in these 


equations gives 


Vir ~ (Vog/VQ)YZTA/M; Mw [ey Ria 


Next, the net energy, dE expended to accelerate (and decelerate ions) to C. 


upon trapping them is determined from 


: : 2 2 
dE f(c hen) (Mic. /2 - Mi (c. - Mace /2) dv... 


2 2 
Bitar el 8 aha ira na ann aay i at 


5 S 


Eq. 4.2 


Integrating the right hand side of this equation gives the total energy ES absorbed by 


the trapped ions. A Maxwellian distribution is assumed for f and the trapping velocity 


View is normalized to View a . With these assumptions, Equation 4.2 
S 
reduces to 
F_ = 


r= “enn (ZT, + 3T;)(ZT/T, + 3) 1/2 , 


fore evE(2 = v) ex(-274/2T, = 3/2)(1 - v)? - 


(2 + V) exp(-2T,/2T, - 3/2)(1 + v)7I} av 


Eq. 4.3 


3 1/4 
vmax (Vo6/Ve) Yu /2y. MT ATT + 37.) R 
The available ion wave energy is a small fraction of the incident electromagnetic wave 
energy. If the backscatter reflectivity grows linearly from O to R, the energy available to 


the ion waves inatime 7 canbe approximated as 


Ey = IRw.t/2Lw, Eq. 4.4 
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Equating EW to the energy delivered to the trapped ions allows the risetime T to 
saturation to be estimated. Figure 4.18 shows the results of such risetime calculations 
for various experimental parameters of interest. What is illuminating is that if steady state 
reflectivity R increases, the risetime does not diminish proportionately because even 
though more energy goes into ion waves, more ions are trapped and heated. Also, it is 
clear that lower density plasmas have fewer available ions to trap and heat and the 
risetime decreases. In the case of the short pulse CO, laser pump where lower electron 
densities were measured interferometrically, faster instability risetimes were detected. 
Even with this crude model, the estimated time scales are close to those observed 
experimentally and qualitiatively substantiate the differences between injection mode 
locked and gain switched CO, laser pumps. Incidently, as the instability growth rate is 

« MEK , one would expect the high density plasma with the gain switched 
CO, laser to exhibit faster risetimes. The opposite trend gives further credence to the ion 
trapping and heating picture presented here, at least until ion wave convection occurs. 


The characteristic time scale for this convection would be = L/c. , which 


‘conv 
is much longer than the instability risetimes. 


4.7.5 Wavelength resolved Thomson Scattering 

The final Thomson scattering experiment was undertaken to measure the 
frequency spectrum of the scattered light. Unfortunately , time broadening through the 
monochromator made high speed streaks unuseable. Instead, time integrated scattering 
measurements were made with the PARC Optical Multichannel Analyser. Absolute 
wavelength calibration was difficult to achieve because of the small wavelength shifts 
from do which were present, although the instrumental resolution ( 0.4 A) was 
sufficient to observe the scattered line width after deconvolution . Figure 4.19 shows a 
typical Thomson scattered light spectrum when the gain switched CO, laser irradiated the 
oxygen plasma (| = 5 x 10}? W/cm’). For comparison, the 0.30 A spectral FWHM can 
be correlated with the infrared backscattered light spectral width. From Thomson 
scattering theory, dr JA, = dv/c * 0 where dd. is the 
scattered ruby light spectral width and dy is the spread in the ion wave phase 


velocity. The Brillouin ion mode wave number is 2k, , so the frequency spread in 
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Figure 4.18 SBS Risetime Estimates from lon Trapping 
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the infrared backscattered light is 2ko dv and the corresponding wavelength 
spread is dv = (2.,/0) : (da /r.) . Thus for the 0.30 A Thomson 
scattered light spectral width, one expects a 70 i infrared backscattered spectral 
width. This estimated dd is larger than the directly observed 30 A average 
backscattered light spectral width. It is likely that di was overestimated since the 
Thomson scattering detection instrumental width was larger than dn, In any case, 


the Thomson scattered light spectral width again indicates strong ion wave damping. 


4.8 Termination of SBS in the Laser-Plasma Interaction 

The fundamental question still remains as to what caused the stimulated Brillouin 
scattering to terminate well ahead of both the IML and gain switched laser pump. 
Interferometry showed that the plasma was still sufficiently dense to provide a high gain 
medium and no anomalous density structure such as profile steepening was seen. Pump 
cutoff through critical layer reflection was also unlikely since the Brillouin backscattered 
light ceased even with O Torr helium background gas where underdense plasma was 
assured. The remaining possibilities include filamentation with resultant electron heating, 
increased pump absorption through enhanced ion turbulence or inverse bremsstrahlung 
absorption, large velocity gradients resulting in wave phase mismatch or plasma 
susceptibility changes. 

Filamentation occurs when transverse pump inhomogeneities intensify and the 
pump wave is channelled into a number of high intensity filaments. This happens because 
the incident light beam’s ponderomotive force on the plasma_ electrons 


eter -(e°/2m.w°)v<E*> , drives the electrons from the region of 


highest ns intensity. The resulting refractive index change focusses the laser light to 
create even higher local pump intensities to reinforce the original pump inhomogeneity. 
Although the light intensity is higher due to filamentation, the decreased plasma density 
and elevated electron temperature can result in an overall reduction of the Brillouin gain. 
When ZT O/T. 2) , the convective Brillouin growth rate is much larger 
than that for filamentation. However, from the results discussed in this chapter, the ions 


can be significantly heated to Landau damp the ion modes associated with Brillouin 


scattering (this is often called Compton scattering). In this case filamentation can compete 
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with Brillouin scattering. The maximum filamentation growth occurs for k, = Q and 


ie 


Bins 2 
lk, = (n/4n dvi civ a) (1/(1 + T,/ZT,)) 
where k,_ is the filament wavenumber parallel to ko and k , is the perpendicular 


component [12]. The spatial growth rate of the filament for large damping is 


Im k, = v1 - n/n. (199° /845C) (Vog/Ve) (1 rere] 7) 


Since the density fluctuation wavenumbers associated with filamentation are not normal 
modes, these growth rates are not diminished by increased Landau damping as they are 
for the Brillouin instability. Brillouin scattering has a spatial growth rate larger by 
OY (w./2y,) which means that for heavy damping, the growth rates of the two 
instabilites are Comparable. For n/n. SHOrbuy, (Vie/Ve)’ = 2 T. = ZT, _ the 
optimum filament transverse dimension is 30 microns which indicates that if a single 
filament forms (whole beam self focussing), the local pump intensity can increase tenfold 
over the usual focal spot intensities. Density rarefaction in the filament through the 
ponderomotive force will act to decrease the Brillouin gain and also the filament 
transverse dimension as calculated above, will enlarge. Further, filamentation can give 
electron heating as well [12], thus reducing (v</vQ)° Although difficult to 
confirm, there is separate experimental evidence for filamentation in the jet (Of fenberger 
[53]) and it is feasible that filamentation could provide a mechanism to inhibit Brillouin 
scattering, at least until normal plasma rarefaction could do so. 
A second possible mechanism to interfere with the Brillouin instability is enhanced 
short wavelength ion turbulence. An increase in the plasma resistivity due to strong 
density fluctuations, results in increased light absorption over classical inverse 


bremsstrahlung (Dawson and Oberman [63]). Faeh! and Kruer [64] show that the ion 


turbulence absorption length is 


toe 4vq/[w,(an/n,)°] = 2,/0¥1 ~ n/n. xfon/n, [7] 
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2 . 
where | n/n | is the square of the normalized homogeneous and isotropic 


density fluctuations. Ruby Thomson scattering measurements on the gas jet target have 


revealed short wavelength ( KAn = 0.5 ) ion fluctuations in the plasma which 

grew as the Brillouin instability decayed (Al-Shiraida et al. [65]). A 240 micron absorption 
é 

length is calculated from the observed | sn/n.,| = 0.04 . Although this length 


is larger than the 80 micron Brillouin interaction region, the integrated effect of the 
enhanced absorption over its extent in the plasma ahead of the focus could be sufficient 
to interrupt the Brillouin instablity at later time as the turbulence develops. The classical 
absorption length for comparison, was =1.7 mm which is even larger still and is 
therefore unlikely to contribute to decrementing the laser pump. 

Although several mechanisms are postulated for ion turbulence generation, one 
predicts an even worse situation if the Brillouin instability itself is participating in the ion 
turbulence generation. lon trapping in the Brillouin ion fluctuation potential can excite a 
bump-in-tail instability (Davidson, [28]) whereby short wavelength ion fluctuations are 
generated (Kartunnen [62]). This new instability is an energy sink to the Brillouin ion 
fluctuations which lowers the Brillouin growth rate and at the same time is a source for 
the ion turbulence. The effectiveness of the mechanism has not been studied thoroughly 
but qualitatively suggests that the short wavelength ion turbulence produced in this 
manner can have a twofold influence on the Brillouin instability through laser pump 
absorption and ion wave damping. The ion fluctuations produced by this instability are 
nominally parallel with La and do not affect absorption. However, mode coupling and 
diffusion in k space may be important for large fluctuations. This bump in tail instability 
also succeeds in explaining the relative timing of the the two instabilities. 

In conclusion, the abrupt termination of the stimulated Brillouin scattering is as yet 
unexplained though it is plausible that filamentation and strong ion turbulence which are 


known to exist, may account for the observed behaviour. 
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5. Conclusions 

A comprehensive study of short and long pulse CO, laser induced Brillouin 
scattering at intensities I ~ O(10%3 W cm”) has been done. The careful characterization 
of the plasma parameters which included electron temperature, electron density and 
electron density gradients, allowed a quantitative interpretation of the stimulated Brillouin 
scattering results. Both the electromagnetic and electrostatic waves of the Brillouin decay 
triplet were detected and the dominant laser and plasma effects were identified. 

The broadened, red-shifted spectrum of the CO, laser light backscattered from 
the oxygen plasma positively identified stimulated Brillouin scattering as the generating 
mechanism. Backscatter reflectivity and spectra for both the gain switched and IML laser 
suggested that a self consistent ion trapping/heating effect was responsible for elevating 
the ion temperature as required for ion Landau damping to broaden the instability’s 
bandwidth. For the gain switched laser, a large ion wave damping coefficient, 

Y,/W, =0.6 was estimated from the infrared backscatter reflectivity. The 
comparison of a kinetic, frequency dependent gain calculation to the backscattered light 
spectra yielded electron-to—ion temperature ratio estimates, ZT /T. =1. This small 
temperature ratio confirmed that strong ion Landau damping took place. 

The IML laser backscatter results were a convincing demonstration of ion heating 
from the Brillouin interaction. The interaction time was too short for significant 
electron-ion collisional equilibration and yet the broad backscatter spectra again 
substantiated an ion trapping/heating description. A simple model prediction based on ion 
heating accounted for the IML laser reflectivity including the high backscatter at low laser 
intensities. The ion trapping/heating model was prepared from the WKB coupled mode 
equations with ion Landau damping from a fraction of the ions heated to temperatures 

0(ZT,) 

lon fluctuations driven by the Brillouin instability were diagnosed through Thomson 
scattering of a ruby laser probe beam. In Thomson scattering wavenumber 
measurements, the detected ion fluctuations were clearly those of the Brillouin instability 
since the measured k wavenumber corresponded to that expected from 
Brillouin wavenumber matching. No Ko or 4k, ion wave harmonics were 


detected in the Thomson scattering experiments, thus eliminating ion wave decay as a 
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Brillouin saturation mechanism at our laser intensities and for our plasma conditions. In a 
second experiment, measurements of the normalized ion fluctuation level, én/n ; 
and its VR reflectivity dependence, verified that convective saturation of the 
Brillouin instability through ion heating had indeed taken place. A third experiment was the 
direct observation of the Brillouin interaction zone through spatially and temporally 
resolved ruby Thomson scattering. The observed 80 micron interaction length was 
explained through wavenumber mismatch of the participating waves away from the 
instability's resonant point. This length coincided with earlier estimates made from the 
infrared backscattered ligh reflectivity. In these same experiments, the 500 ps ion 
fluctuation risetime was strong evidence for ion heating. These slow risetimes agreed 
with the model of heated ions diverting energy from the Brillouin ion waves to lower the 
Brillouin growth rate. Risetime estimates from the ion trapping/heating model were close 
to those observed. 

A successful attempt was made to reduce the Brillouin scattering by using a 
multiwavelength CO, laser. The selected laser lines had to have frequency separations 
larger than the Brillouin homogeneous growth rate otherwise anomalous mode coupling 
occurred and reduced the effectiveness of the multiwavelength laser. This unique 
demonstration of SBS reduction is encouraging evidence that simple laser modifications 
can potentially improve laser—plasma coupling in fusion experiments. 

The termination of SBS in the laser-plasma experiment is still not fully 
understood, although filamentation and enhanced pump absorption remain as the most 
likely candidates to lower the SBS gain below threshold. Both filamentation and strong ion 
turbulence have been observed in other experiments made with the CO, laser and gas jet 
target. Indeed, the cessation of SBS just before enhanced ion turbulence takes place, 
suggests that strong laser absorption could be interrupting the Brillouin instability. 

In summary, the complete characterization of the laser, plasma and interaction 
parameters permitted the quantitative interpretation of the observed stimulated Brillouin 
scattering. Spectral, spatial and temporal measurements of SBS and the associated ion 
fluctuations strongly support the model of ion heating and wave damping along with the 


finite matching length. 
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Appendix A 


Nonlinear Equations 
Several nonlinear equations had to be solved including the ion heating reflectivity 
calculation in Eq. 2.26, the ion trapping in Eq. 2.32 and the stimulated Brillouin dispersion 
relation in Eq 2.11. The first two nonlinear equations were solved through Newton 
Raphson methods implemented on a computer. Except in the simplest nonlinear equations, 
the function derivative was approximated by a forward difference so the iteration 


formula became: 
=x, - hf(x,)/(f(x, + h) - f(x;)) 


Careful selection of initial values of X were required due to the exponential 
behaviour of the nonlinear equations. The simultaneous nonlinear equations describing the 
stimulated Brillouin dispersion were solved through an iterative substitution method as 


described in Chapter 2. 


Ordinary Differential Equations 
Systems of ordinary differential equations had to be solved for the rate equation 
model of the injection mode locked laser, the spatial coupled mode equations, and the 
temporal coupled mode equations. A second order Runge-Kutta integration formula was 
used with interval halving to ensure that the error was within 0.5%. The two step 
computational method used to calculate y; ee oc IN was: 


= Vii + hy" (X5 Vy gore yy )/2 
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The above shows the calculation of Vere from as with step size 
riaet sag 
h 


Deconvolution of Spectra 
The analysis of spectral data was complicated by the instrumental function of the 


dispersing elements. In particular, the resolution of a grating spectrograph is determined 


by the number of grating lines, N , dispersing the incident light, ie., 
A\/rA = 1/N for first order diffraction. For infrared light measurements 
with the 75 lines/mm grating spectrometer, the resolution is 28 A which is the 


same order of magnitude as the detected line shifts and width in Brillouin backscattered 
light. The pseudodeconvolution of Jones et al [47] was used to remove the instrumental 
responses from the spectra. This iterative method involves convolving the original data 
with the instrument response and the resulting function is compared with the original data 
to give a measure of the spectral broadening. The algorithm for this discrete point 


method is: 


(1) S' = S. Xx (S./C(S.)) 


POV Ges ten 
= S' x (S,/C(S')) 


where C(S;) is the convolution of the j th spectrum ce 
J 
with the instrument function. To reduce computational noise, a five point least squares 


smoothing was passed over each S. 


Sok oie ikke) ee ms on (kezyen akki2) 2 


where S; k is the k th averaged point of S. Typically four 
iterations were used to deconvolve the spectra and Sy were convolved with 
the instrument function for comparison with the original data So to ensure that 


it could be recovered. 
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‘ Abel Transform 
The computation of radial functions from data integrated along cords in a 
cylindrically symmetric system is achieved through the Abel transform. This transform 
was applied to a numerical routine used to invert plasma interferograms and obtain the 
radial electron density profiles at various axial positions in the plasma. The usual inversion 


formula for refractive index can be written 


=1 {U's "ds -1/2 
u(r,z) - 1 =— — (u - U) du 
™ Jy du 

where u(r ,z) is the refractive index at radial position Yr in a cylindrically 
symmetric plasma of radius ce The fringe shift at u is § where 
u = r2 and Uo = eee . The refractive index profile is related to the 
electron densit b ae eae Z 2 from which 

y y u noe fe Moro by 
na (tz) is found. A linear approximation to ds/du between data 


points reduces the integral to the numerical computation: 


N-1 Sowers 
=, 3 j+] j 1 
a | a 2 ~——4 [(y +] Sen) a (u. u.)'/25 
j=l ea Re 50 es J 1 J i 
bes 
The refractive index at the j th of N equispaced data points is 
ue . Although this is a highly simplified computation scheme (for comparison 


see for example i66]), tests done on analytic solutions confirmed its accuracy to within a 


satisfactory 5%. 


Plasma Susceptibility Calculation 


The kinetic limit of the SBS dispersion relation in Equation 2.11 required the ion 


susceptibility X Maer S (ky 2/k?)W(z) where 


Z 
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and z= w/V-k which is complex. The function W(z) is 


related to the usual plasma dispersion function [(z) (Fried and Conte [17]) by: 


W(z) =(z/72)+r(z/v2) +1 


Two computations were appropriate for calculating W(z) and T(z) 

depending upon —_Im(z) . When |Im(z)| > 1 , a continued 

fraction algorithm (Fried and Conte) was used to compute T(z) from which 
W(z) was derived. A convergent series (Ichimaru [18]) gave W(z) directly 


for all other cases. These two algorithms are summarized below with corrections to the 


original references. 
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Appendix B 


Free Jet Program 
The method of solution for free jet flow is well known from hodograph theory 
(Shapiro [45]). Reference should be made to elementary texts on supersonic flow theory 
to assist in understanding the equations in this appendix. Figure B.1 is a diagram of the 
free jet showing the Mach lines a-h separating flow fields 1 to 8. Upon crossing a Mach 
line, the gas flow parameters change, consistent with the various isentropic relations. The 


free boundary angle, on is found first by calculating the Mach number , M in 


5 
region 3 from pressure balance Pe = Pa where Po is the pressure in 


region 3 and P is the ambient pressure (pressure difference does not occur across 


A 
stationary free boundaries) and then applying the Prandti-Meyer equation. The change of 


the flow. direction when going from region 1 to region 3 is 


oF = Wy, ~ Ws . The two equations for calculating the above are 
M = L(2r(y - 1))(PyfPy VOY ay ¥? Eq 6.1 
”, yer lil/2 y-l 9 bie 
w(M) = Atan(M™ - 1) - Atan (Mo - 1) 
Ne yt 


EG eZ 

Equation B.1 is an isentropic flow relation and B.2 is the Prandti-Meyer equation. The 
ratio of specific heats, y is equal to 1.4 for a diatomic molecule. For example, if 
PS =2000: Torr, P, =10 Torr and M, =2.96 (from the nozzle calculations 


A 


outlined in Chapter 3) then M. = 4,2] and 


ners -(w(4.21) - w(2.96)) = 68.45 - 48.98 = 19.47° 


Depending upon whether a Mach line is a compression (Type |) wave or 


rarefaction wave (Type Il), the flow direction after crossing them is determined from 
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Figure B.1 Flow fields and Mach Lines in a Gas Jet 
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ts) and ts) defined as 


6, = -w(M) + 2I - 1000 


I 
rr = w(M) + 2II - 1000 
The integration constants 2A= 7 L000 and <¢lI = 1000 


are useful during computations. From above, if one travels from flow field i to flow field 


j, then 

FF Ce S Ws = W, for Type I 

Seg -(w, - w,) for Type II pais 
A type Il Mach line is inclined to the jet axis at an angle Ce ae 
where 


a = Atan(1//M2 - 1) 


Clearly a _ is the angle that the Mach line makes to the downstream streamlines and 


8 is added to it to restore it to an angle with respect to the jet axis. Similary, a type 


| Mach line is at angle oa to the jet axis. Equation B.3 can be rewritten as 


6 = I+ II - 1000 
w= I1- Il 
The variable | is unchanged when crossing a type Il wave, and Il is unchanged when a type 


| wave is crossed. It is then a simple matter to calculate 96 Beet for different flow 


regions in the jet From w , M is calculated with Equation B.2 and then & 
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from above. Once all 6 and a are known, then all Mach line angles are 
determined. From this information, trigonometric relations are used to locate the 
intersection of the Mach lines and the flow field pattern can be drawn. 

The free jet program below is written for a Texas Instruments 99/4A computer. 


The key variables are: 


A(2,8) = (Ill) values for the 8 fields 
OM(8)= w 

TH(8) = 0 for the 8 fields 
MM(8) =  M for the 8 fields 
AL(8) = fe) for the 8 fields 


A1-G1 = Mach line angles to the jet axis 
X X(N), YY(N) = x,y coordinates of the points L to Q 


100 DIM A(28) 

110 REM FIRST CALCULATE W(M1),W(M3) 
120 INPUT "PO,PA": PO,PA 

130 M3=SQR(5*((PO/PA)!(.4/1.4)-1)) 
140 M4=M3 

150 REM 1170 IS W(M) SUBROUTINE 
160  GOSUB 1170 

170 Wi=W 

180 M4=2.96 

190 GOSUB 1170 

200 REM INITIAL AND B.C. 

210 OM(1)=W 

220 TH(1)=0 

230 TH(3)=W-W1 

240 TH(2)=TH(3)/2 

250 TH(4)=0 

260 TH(6)=0 

270 TH(7)=0 

280 REM START OF FIELD CALC. 
290 A(1i,1)=500+0M(1) /2+TH(1) /2 
300 Teo ence ae eee 
31.0 we ildene aA Ging W) 

320 OM(2)=2*A(1,1)-TH(2)-1000 
330 A(2,2)=(1000-0M(2)+TH(2)) /2 
340 A(1,3)=A(1,2) 

350 om (3)= 2*A(1, 2) -TH(3)-1000 
360 a(2,3)=(1000-0M(3)+TH(3)) /2 
370 A(2,4)=A(2,2) 


ee asinns?ac jos, its tate 
ori? “Sissoh bt ‘baal aie, te 


rte seed: via 
rsiuerrios AEBS =InSmBsbzl east ‘oF ee 


; <a 


OM(4)=1000-2*A(2,2) 

ACLs 6h 1 000+om (244TH 4)) /2 
A(1,5)=A(1,4) 
A(2,5)=A(2,3) 

om (5)= AGL, eat, 3) 
TH(5)=A(1,4)+A(2, 3)-1000 
oM(6)= 1000+TH(6) -2*A(2, 5) 
A(2, en A(2,5) 

A(1,6 =1000+0M(6)+TH(6)) /2 
A(1,7)=A(1,5) 

om (7)= 2*A(1,5)- -1000 
A(2,7)=(1000-0M(7))/2 
A(1,8)=A(1,6) 
A(2,8)=A(2,7) 

oM(B)=A(1, 6)-A(2,7) 
TH(8)=A(1, eaters! 
FOR I=1 TO 8 


MM(I)=M 
sur Je =57. 2957 5eAtH(1/SQn(we 1)) 


eae "MACH + STREAM ANGLES" 
FOR I=1 TO 8 

PRINT OM(I),TH(I) 

NEXT I 

PRINT "MACH NO, ALPHA" 

FOR I=1 TO 8 

PRINT MM(I),AL(TI) 

NEXT I 
A1=(TH(1)+TH(2)+AL(1)+AL(2) 
B1=(TH(2)+TH(3)+AL(2 ae 
C1=(AL(2)+AL(4) -TH(2)-TH(4) 
Di= 0 3Yoan3)-m 3-3 
E1=(AL(3)+AL(5)-TH(3)-TH(5) 
F1=(AL(5)+AL(6)-TH( 5)-TH(6) 
Gi= Perera a cats 
Hi=(AL(7)+AL(8)-TH(7)-TH(8 
Q=1/57.29578 
sen 

YY(1)=1 


) /2 
y 
y 
)/2 
) /2 
)/2 


yy (3)=1-(XX(3)-XX(1)*TAN(Q*C1 
XX(2)=XX(3)+(1-YY(3)) /TAN(Q*D1) 
YY(2)=1 
=-XX (4 


yy(4 *TAN( -TH(3)*Q 


Se apie 


te Se il ite nen USES ee 


XX(5)=(1-Y¥(4)+XX(2)*TAN( Q*F1)+XX (4) *TAN(Q*G1) / 


(TAN(Q*F1)+TAN(Q*G1 )) 
ae A eed re ea vn as 
XX(6)=Xx(5 

YY(6) =Y¥ (4) 

PRINT "CHARACTERISTIC COORD." 


+(Y¥(5)-Y¥(4) ) /TAN(Q*H1 
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900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 


1060 
1070 
1080 
1090 
1100 
4210 
1120 
1130 
1140 
43:50 
1160 
1170 
1180 
1190 
1200 
1220 
1220 
1230 
1240 


FOR I=1 TO 6 

PRINT XX(I),YY(I) 

NEXT I 

PRINT "FIELD DENSITIES" 
FOR I=1 TO 8 


G1.4 
GG=SQR(G+1) /(G-1)) 
M=1.1 

DM=.0001 

M1=M 

GOSUB 1130 


F1i=F 
M1=M+DM 
GOSUB 1130 
M2=M 
=M2-DM*F1 /(F-F1) 
IF ABS(M2-M)>1E-4 THEN 1040 
RETURN 
MR=SQR(M1*M1-1) 
F=57.29578*(GG*ATN(MR/GG )-ATN(MR) )-T 
RETURN 
REM W(M) SUBROUTINE 
GG=SQR((2.4/.4)) 
MR=SQR(M4*M4-1 
W=57.29578*( (GG*ATN(MR/GG) -ATN(MR) ) 
RETURN 
G=1.4 
P=PO/(1+.2*M*M)!2.5 
NC=P/2000*12*6 .41 
RETURN 
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Appendix C 


Fresnel Mirrors Interferometer 

The plasma electron density measurements reported in this thesis were generated 
from the Abel inversion of interferograms. Implementing a workable and convenient 
interferometer was not as trivial as it would first seem. The usual folded wavefront and 
Mach-Zender interferometers were tried and found unsatisfactory. Both of these are 
division of wavefront interferometers with an unacceptable 50% reduction of the 
interferogram exposure. Also, the Mach-Zender optical elements were located inside the 
experiment target chamber, making adjustments awkward. Finally, the minimum number of 
critical components in the Fresnel Mirrors interferometer made it attractive to use. 

Figure C. 1(a) and C.1(b) will be referred to in explaining the interferometer. Several 
design goals had to be Satisfied simultaneously through the appropriate selection of 
lenses’ focal lengths and positions. These design objectives included the image 
magnification, fringe spacing, lens L1 clearance from the CO, laser beam and not 
damaging the exit etalon with focussed Ruby laser light. The expanded and collimated 
laser light passed through the plasma position O1 in figure C.1(a). One half of the laser 


light was the reference beam and the other half the test beam. The telescope lenses were 


separated by f, + fh to recollimate the Ruby light. At the same time the 
fringes were localized at the plasma. From the thin lens equation, the rea! image location, 
I, and telescope magnification M were 
* * 

= = 2 Fo: Gal 

I, (1 +M )fo 01M g 

M = I,/0 

ae Eq. C.2 


where M* = Fo/F, 


The overall magnification, M was the product of M and the camera 


t 


magnification, Mion . The Fresnel mirrors were slightly inclined to each other 


(angle ¢ in Fig. C.1(b)) to cause beam overlap. By locating the mirrors very close to 


LZ: Dw I, . Further, the interference fringe spacing at I, was 
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Figure C.1 Principal interferometer Components 
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* 
AX = AT, /M , where A was the Ruby laser wavelength in this case. 
* 
The width of the overlapping area was yM . Referring this fringe spacing back 
to the plasma, the interferometer’'s line resolution is then: 
AX. = Ax/M 
p Eq. Gc 
Some convenience results in normalizing O1 and obtaining slightly different expressions 


for the magnification and resolution. These are: 


N= iM + m2(1/0,- Elia 


Eq. C.4 
ax. = 0, f / re 
sh leans eg BA Eq C5 
* 
h = 
where 0, 0,/f, ‘ z 
Figure C.2 shows 0, as a function of M for various M 
* * 
Once M and M_= are selected, 0, is easily determined from Fig. C.2. 
* * * 
Also, I, = MO, where I, = I,/f, . One set of parameters for the 


interferometry are indicated with intersecting dashed lines. In this particular arrangement, 

f, = 30 cm and rm = 0.7 cm, and the resolution is calculated with 
Equation €.5 to be 26 microns. This fringe spacing corresponds very closely to the 25 
micron resolution estimated from counting 8 He-Ne laser fringes across a 200 micron 
diameter alignment microball positioned at the CO, laser focus. The overall magnification, 

Mo = Mx We ey was 50X onto the Polaroid type 410 film in the 
interferometry camera. It should be mentioned that CO, laser beam clearance and a safe 


Ruby laser intensity onto the exit etalon were also satisfied in this configuration. Normal 


Abel inversion analysis was then done on the collected interferograms. 
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Figure C.2 Interferometer Design Parameters 


; 


ie 
iy " 


ul 


| ie ; 
é my pA ay 


i) 
haan 


oF lt J a 
Cah ay, 
Pl ke 

Cm) Poy ql 
rah Mees, 


